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ABSTRACT 


The equations for estimating the performance of an axial flow 
turbine are examined to determine the efficiency and other parameters 
when very high power outputs per stage are required. It is indicated 
that turbine stages should be designed for peak efficiency or for 
high power output, but that compromise designs are not advisable. 

A survey of the immediately available literature on losses in 
turbine stages at the U. S. Naval Postgraduate School is outlined. 

An attempt to find and correlate specific loss data with theoretical 
relations was not successful. Data is fragmentary, theoretical re¬ 
lations are few, and the correlation of the available data on profile 
losses, tip clearance losses, effect of trailing edge thickness and 
secondary losses is poor. 

An axial flow single stage turbine was designed for transonic 
flows and a digital computer program was written to evaluate the 
off-design performance of this turbine. The particular relations for 
supersonic expansion after a blade row are included in the program 
as well as the effects of entry shocks into the rotor. Initial 
indications are that the program should be useful for evaluating 
turbines designed by the procedure shown. 

The computer programs included are written in FORTRAN language 
for the CDC 1604 Digital Computer. 
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THE PERFORMANCE OF HIGH POWER OUTPUT AXIAL 


FLOW TURBINES UTILIZING FLOWS IN THE TRANSONIC REGIME 

I. Introduction 

The increasing need in space applications for turbines of high 
power output from a minimum physical volume and weight requires more 
knowledge of efficiency, losses and performance estimation of these 
turbines. 

If the required power output and space limitations are predominant 
then efficiency must suffer. As a means of showing the degree of 
efficiency sacrificed, the equations which determine the power output 
and efficiency of an arbitrary turbine were programmed for digital 
computer solution in terms of the head coefficient. The head coef¬ 
ficient for peak efficiency was determined and then the head coef¬ 
ficient was increased in increments and various parameters affecting 
the turbine performance were computed at each increment until the 
efficiency dropped to an arbitrary cutoff point. The data are plotted 
in non-dimensional form and simple equations are found to represent 
the data. 

Of paramount importance in the design of a turbine is the esti¬ 
mation of the losses due to boundary layer formation, mixing effects, 
clearances and secondary flows. This area is explored in the hope 
of finding a more accurate and precise method of presenting these 
losses. A survey of the immediately available literature is included. 


A particular single-stage turbine is designed for supersonic 
flow leaving the stator blade row and a relative subsonic flow as 
seen by the rotor. The equations for determining the effects of 
after-expansion behind a blade row and entry shocks entering a 
blade row are programmed for the computer in addition to the 
regular performance parameters. The rpm and pressure ratio across 
the turbine are then varied to show the effects of demanding extreme 
power outputs from a turbine so designed. A performance map of the 
turbine is developed. 


-2- 


II. Development of High Head Coefficient Data 


A. Basic Equations 

In Ref. 1 the equations are developed to establish the design 
parameters for turbine stages. These equations are obtained by 
assuming that the flow through the turbine is represented by the 
conditions on a mean stream surface and the primary parameter is 
taken to be the "head coefficient". The head coefficient relates 
the total isentropic drop of the stage from the inlet stagnation 
pressure to the discharge static pressure to the kinetic "energy" 
of the peripheral speed at the rotor entrance and is defined as 



(i) 


In order to relate the theoretical expressions to values found 
in engineering use, some expression for the losses involved must be 
used. A velocity coefficient (^) is used as a measure of rotor 
efficiency and is defined as 



( 2 ) 


This velocity coefficient has been found through experience to depend 
primarily on the deflection of the flow in the rotor (/\A ) . Vavra 
has established a mean curve through data from a variety of sources 
to be 




4 ^ 97 _ 

/80-A/S 


(3) 
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and this equation is used to represent the losses for purposes of 
evaluating the efficiency at each step of the calculation. 


Following the procedure outlined in Ref. 1, arbitrary values 
were chosen for the angle of the flow leaving the stator (Qf^), 
reheat factor (f), and the carryover coefficient of kinetic energy 
into the next stage , the radius ratio (R 2 /RP , the meridional 

velocity ratio (V j/V j), the velocity coefficient in the stator if). 


and the carryover coefficient of kinetic energy from the stator into 

V 


the rotor ((p ). Then the degree of reaction (r ) is selected and 


the head coefficient is increased in increments. The degree of 
reaction is the fraction of the isentropic enthalpy change through 
the turbine stage which occurs in the rotor. At each increment of 
head coefficient the particular values of velocity coefficient , 

internal efficiency (^ ), leaving loss coefficient (k^), ratio of 
meridional absolute velocity to peripheral speed (V^/U^) , relative 
velocity ratio (W^/W^), and the relative and absolute angles (0^ > 

^ 1 ? jS 2 ) are computed. 


A similar computer program using the same equations was written 
to find the peak efficiency for a given entering angle and degree of 
reaction, and then all the other parameters are computed. The values 
are also found for higher head coefficients such that the efficiency 
is 0.5% less than optimum and then 1.07 o less than optimum. 
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B. 


Computer Program 

The equations to be solved by iteration until *4 are; 


1 . i ■■ - 


— — 


/ 


f £<XLO(/ i7-A* 




(A 


AT 


where 


(Vmz/^n/) /U&o(/ ''J/—/xj^ 

>2 I 

-2& 


(A) 


(5) 


(5a) 


3 


4.97 


/O* r /80 A /4 


(3) 


Once this condition ( ^ 3 ^ ) is satisfied the values of the 

A 

other parameters are found from 


1 . 


* - m[ f ^ ^wM + Hf) 

-J 

4 = -4 t z x z 


( 6 ) 


(7) 



'Sl 


{dZ "'U/ + ^IK 


/ 


2 1 


c 

u, 


= y //T? /£/ Coart, 


( 8 ) 


-5- 






































-A _ _ (&/&)/&£ ^^ 

A - 


(9) 



The equations were coded in Fortran language, Ref. 2, and solved 
on the Control Data Corporation 1604 computer at the U. S. Naval 
Postgraduate School. A copy of the program is presented in Appendix III. 
The resulting data printouts are included in Appendix V as Table I. 


The program finding the peak efficiencies, efficiencies 0.5% less 
than optimum and efficiencies 1.07o less than optimum is also in 
Appendix III, with data in Appendix V as Table II. It should be 
noted here that all equations are completely coded in Fortran language 
so that any selected numerical values may be introduced in the be¬ 
ginning of the program. Any of these arbitrarily chosen values may 
be changed and the equations of the main program remain unchanged. 


C. Results 

Representative data from Tables I and II are presented in Fig. 1. 
The regularity and progression of the data curves indicate that non- 
dimensionalizing the efficiency and head coefficient parameters might 
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Fa profitable. When all the efficiencies for one entering angle and 

one degree of reaction are divided by the optimum efficiency for that 

set of conditions (/^/ / 0 ^^} an< ^ the head coefficients are divided 

by the optimum head coefficient (k. /k. ) the curves of Fies. 2 . 

is is opt ° ’ 

3, 4 and 5 result. 


A closer inspection of these plots reveals that for each degree 
of reaction, one curve suffices for all four entering angles. 

The cubic equation 



-.000223 


l k is ^ 

3 

+ .008345 1 

I } 

' is 

2 

1 O A/ O 1 

f k i s \ 

k. 

^ is opt 

k 

[ is opt 

-.12068 

k. 

1 is opt 1 


+1.125 


represents all the data for r “ 0.0 for all entering angles from 60° 
to 75°. 


Similarly as an approximation to all the data at r = 0.25, the 
equation is 


<?/? opt > * " • 000223 


k. 


is 


is opt 


+.00770 


k. 

is 


k. 


1098 


k, 

is 


is opt 


k. 

is opt 


+ 1.1 


and for r * 0.50 



-.000125 


/ k. 
is 


k. 

is opt 


3 

+.00554 


is 


k. 
is opt 


2 / k - 

-.09414 -- 

1 is opt 


+1.097 


Ihe maximum error in the range of entering angles presented here is 
approximately 8 parts in 130. This is probably accurate enough for 
preliminary engineering estimates based on the equations in Ref. 1 
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used in the computer programs. As will be seen later, the loss data 
is not known to this accuracy. 

The curves show that, roughly speaking, increasing the head 
coefficient to five times the optimum causes a reduction of 307> of 
the obtainable peak efficiency. Increasing the head coefficient to 
ten times the optimum causes a reduction of 457o of the obtainable peak 
efficiency. This indicates that if it is necessary to go to very high 
head coefficients to obtain a certain required work output, the rate 
of change of efficiency with increasing head coefficient is greatest 
from optimum to five times optimum head coefficient. The rate of 
change of efficiency decreases as the head coefficient is further 
increased. 

Table I of Appendix V can be used to establish the other stage 
parameters after the head coefficient range of interest is chosen. 

Figs. 1 through 5 are not the off-design performance of a 
particular turbine, but they are the locus of points of an infinite 
number of different turbine designs. Each point on each curve 
represents a turbine designed for that condition. Therefore, if the 
requirement is that the turbine be designed for a high head coef¬ 
ficient, the loss in efficiency need not be the absolute controlling 
factor, and the design may as well be based on extreme high values 
of head coefficient. These investigations show that it is advisable 
to either design the turbine for peak efficiency of maximum output 
but do not compromise both by operating at only two or three times 
the optimum head coefficient. 
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III. Loss Investigation 


Since the design of a turbine depends on the loss assumptions 
made, the use of an average velocity coefficient based on general 
considerations of flow deflection in the rotor is not a very precise 
method of estimating the losses. This was the method used in the 
preceding section of this thesis. A better way should be found and 
the following considerations are basic to the problem. 

A. Loss Components 

It seems logical to class the losses encountered in turbines 
into four particular categories; 

1) Profile loss -- the friction loss associated with the 
formation of the boundary layer on the blades of the cascade. 

2) Tip clearance loss “ the pressure and energy loss associated 
with the leakage flow between the tips of the blades and the casing 

of the rotor, including the effects of trailing vortices at the hub 
where there is zero clearance. 

3) Trailing edge loss -- the loss due to mixing effects behind 
the blades of the row as caused by the finite thickness of the blade 
trailing edge. 

4) Secondary loss -- the loss due to the vortices which are set 
up in the flow in the curved channel between blades giving a velocity 
component normal to the desired direction. 
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In the axial flow turbine, the closely packed blades and large 
flow deflections usually require the flow through the turbine to be 
considered from the point of view of "channel flow" between blades 
instead of using airfoil theories. As will be seen in the subsequent 
discussion, this secondary flow causes the principal problem pre¬ 
venting the accurate assessment of the losses. 

Looking at the flow between two blades of an arbitrary cascade, 
Fig. 6, the velocity distribution is as shown for the mid-span of 
the blade channel and it can be seen that a static pressure gradient 
is developed between surfaces 0-1-2 and 0 1 -1*-2 1 . The centrifugal 
forces on the fluid being turned cause an increase in the static 
pressure near the O'-l 1 -2* wall and decrease the static pressure 
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Assuming non-viscous flow outside the boundary layer requires the 
flow to increase velocity near the 0-1-2 convex wall to maintain the 



Conversely, the flow 


total 


velocity must decrease near the O'-l 1 -2' concave wall outside the 
boundary layer. In the boundary layer at the top and bottom of the 
channel, the loss of velocity (and total pressure) due to friction 
causes the flow direction to be dominated by the pressure gradient 
imposed (assuming the static pressure is constant through the boundary 
layer). Therefore, as the flow progresses from station O-O 1 to 
station 1-1* to station 2-2 1 , the fluid in the upper and lower 
boundary layers tends to flow from 0* toward 1 and from 1' toward 2. 

This change of direction of the flow in the boundary layers 
results in the formation of two vortices superimposed on the main 
flow leaving the channel at station 2-2*. In Fig. 7 the flow patterns 
in the upper and lower boundary layers are shown as well as the 
induced vortices rotation direction at the exit plane. 

This vorticity results in velocity components developed perpen¬ 
dicular to the desired exit velocity. Since a velocity component 
normal to the desired direction cannot be recovered usually, most of 
the useful energy in these components is lost in friction. This is 
the so-called "secondary loss". For a rotating row of blades, the 
blade tip clearance also has a marked effect on the secondary flow. 

The leakage from the high pressure to the low pressure side of a blade 
encourages an additional trailing vortex to form a reaction blading. 
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For impulse blading the flow over the tips which is not deflected also 
enhances vortex formation. For cascades of infinite blade height 
(two-dimensional only) there is no secondary flow since the mechanism 
depends principally on the annulus boundary layers. 



With the physical parameters of blade height, trailing edge 
thickness, tip clearance, profile shape, and spacing given, it should 
be possible to make a reasonably exact estimate of losses. 

B. Survey of the Literature 

t 

The immediately available literature at the U. S. Naval Post¬ 
graduate School was surveyed to find the particular loss breakdown 






































































used by the various sources, the theory and representations offered, 
and any precise loss data published. The survey is limited in scope 
and critical since the particular desire is to locate specific loss 
theories and data. The results of this survey are presented here 
in the order of reference books, general texts, and then periodicals 
and reports. The NACA/NASA reports provide a wealth of data on the 
efficiency, or total loss coefficient, of specific turbine designs 
from tests. Very little information was found for loss components, 
however, in the frame of reference considered here. 


1. In Ref. 1 an equation is derived for the determination of 
a flow coefficient in a rotor as a function of the change in 

flow angle of the relative flow. This equation is a fitted curve 
through data from several sources and a variety of blade shapes and 
operating conditions and is the same equation used in the preceding 
section of this thesis. In Fig. 8 is shown the general shape of the 
curve and the derived equation (3). 


A loss coefficient can be defined as 


Considering the case of * 0, we see the loss coefficient 

f 0 - » -t a represents those losses in the cascade which are 
present when the flow is not deflected. As the flow is deflected 
an amount &, the reduction in ^ indicates an additional loss. 
The additional loss will be taken as the secondary loss. The 
secondary loss coefficient is then 




(id 


( 12 ) 
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Curves shown in subsequent discussions of secondary losses will now 
show a loss coefficient curve based on this equation plotted as a 
dashed line and labeled "VAVRA". This curve will act as a general 
reference and will show "average" conditions so a comparison of each 
method can be made. 



Fig. 8 Flow Coefficient ^ as a Function 
of Turning Angle in a Cascade 

2. A veritable fountainhead of theory, design and ideas are 
presented in Ref. 3. This basic work by Stodola touches on almost 
all the subjects considered important now. The important works of 
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others up until 1927 are included to show their measurements and de¬ 
velopments. The losses through fixed rows of blades are given in 
the form of the so-called velocity coefficient P ; the ratio of actual 
velocity to theoretical velocity. The measurements of the Institution 
of Mechanical Engineers of England and others all show good correlation 
for flows up to sonic speed for stators with an axial entering angle 
and an exiting angle of 70°-78° from the axial. For supersonic 
speeds the data show some scatter due to difference in measuring 
apparatus and some basic differences in design of the convergent- 
divergent passage* 

Secondary flows in bends are explained but there is wide vari¬ 
ation in the available data for friction losses in bends. No pub¬ 
lished data were available for the losses in rotating blade rows. 

The data and experiments on fixed rows of impulse rotor blades were 
of limited application and only show general trends. 

For design purposes, velocity coefficients and loss coefficients 

v 

are given in general terms only, i.e., "velocity coefficient for 

i 

stators = 0.975-0.92 for long and short nozzles respectively". The 
curve presented by Stodola for the velocity coefficient for rotating < 
blade rows as a function of turning angle is included in the data 
used for the determination of the equation (3) by Vavra. 

The tip clearance losses had been explored in a limited fashion. 
These losses were accounted for by an empirical formula for the de¬ 
crease in internal work as a function of clearance and blade height, 
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( 13 ) 


^ =• 6.27 k 1,4 /h 

where k * radial clearance and h is the blade height corresponding 
to zero clearance. This relation, credited to Anderhub, is 
plotted on Fig. 9 and noted as . 

The steam turbine design procedures presented are strongly 
dependent on making the new design a small modification of an already 
designed unit which is in service and for which the efficiency 
is known. The peak efficiencies were recognized to be in the regions 
of k^ g from 1.0 to 5.0, as shown in the previous section of this 
paper. 

3. The axial turbine section of Ref. 4 gives a lucid and brief 
explanation of the flows in such machines but the loss data and 
secondary flow considerations are based completely on Ref. 22, 23 
and 24 which are considered later in this paper. 

4. Since books written to be used as school texts of necessity 
try to cover the complete spectrum of fluid mechanics, theory, and 
design, only the most cursory reference is made to specific losses 
in turbine stages. Refs. 5, 6 and 7 are representative of books of 
this type. Ref. 5 defines a work recovery coefficient as the ratio 
of useful work extracted by the rotor to the total kinetic energy 

of the flow leaving the stator, times the stator flow efficiency. 

This work-recovery coefficient variation with the velocity ratio 
(U/V^) is shown for impulse and reaction conditions and no other 
explanation of the nature of the losses is offered. The charts 
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presented are based on data from Ref. 16, in which report the losses 
in a blade row are presented only as a function of stagger angle and 
incidence. 

In Ref. 6 the author makes the statement that no adequate 
criterion is available for the difference in the flow conditions 
between an efficient and an inefficient hydrodynamic machine of a 
given type. Any losses and efficiencies of turbines offered are 
tied to airfoil theory or definitions of losses by Stodola and 
others. Secondary flows and radial flow of the boundary layer are 
mentioned only from Ref. 17 where the measurements were made on 
widely spaced blades on a rotating blade row at low speed. 

A large list of references on the work in gas turbines is 
given in Ref. 7 for the data available at that time. Principally 
the work in Ref. 16 below is mentioned in connection with turbines 
but only general terms are offered and no specific loss derivations 
are offered. 

5. The most comprehensive effort to make a breakdown of losses 
and show specific data found in a textbook is that of Ref. 8. This 
book takes selected data from Refs. 15, 16 and 23, discussed later 
in this paper, and presents plots of the losses due to leakage, 
profile losses, losses from turbulence and wall friction (including 
secondary loss effects) and incidence losses. The discussion is 
brief, however, and no new data or theory are presented. 
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6. An attempt to mathematically tackle the problem of the 
secondary flows in bends is presented in Ref. 9. The equations of 
motion of an ideal fluid are used and frictional effects are ignored. 

An induced drag pressure drop coefficient is defined to account for 
the power loss due to secondary flow. The method is applied to the 
flow in the corner of a 3' x 4' wind tunnel which has two turning 
vanes to simulate a blade row and the method shows the proper trends. 

A quantitative value derived for the loss is computed and found to 

be about twice the measured loss. No other example of the theory 
being applied is available. 

7. A qualitative investigation of the flow in boundary layers 
and wakes of blade rows is offered in Ref. 10. An excellent discussion 
of the general nature of these flows is accompanied by the statement 
that a quantitative theory for these effects is not available. 

8. An analytical method of estimating turbine performance 
described in Ref. 11 is based on losses due to incidence angle and 
losses due to flow of the fluid through the blade row. A blading 
loss parameter is used to represent the losses but only a general 
range of values of the parameter are given. No method is presented 
for determining the magnitude of the blading loss parameter before 
the turbine is built. 

9. The general state of the art to 1948 is described in Ref. 12. 
The paper is not intended to be used as a basis for other than general 
design considerations, however, and some empirical relations presented 
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were attempts to cover both axial turbines and compressors. This 
necessarily limited the worth of such relations for large flow 
deflection turbine loss estimations. Those relations presented for 
strictly turbine use are repeated in Refs. 24 and 25 which will be 
covered in more detail later in this thesis. 

10. The losses in stators are, of course, easier to measure 
than the losses in rotors. Some excellent tests and measurements 
are presented in Refs. 13 and 14. These show that, in general, a 
properly designed nozzle for a turbine will have a velocity coef¬ 
ficient (ratio of actual velocity to isentropic velocity) value of .96. 
This single value is found to be uniform over a considerable range 

of designs and conditions. Off-design values do fluctuate, partic¬ 
ularly in the case of supersonic nozzles but the reasons do not seem 
obscure and loss estimations for stator rows should be on relatively 
firm ground. 

11. Ref. 15 is mentioned as a reference work in many textbooks 
on turbomachinery. This compilation of the then current design 
practices only shows the losses as a function of axial length of the 
blades of the rotor, turning angle in the rotor and incidence angle. 

The curves for loss as a function of turning angle are included in 
the data used for the determination of the equation (3) by Vavra. 
Specific detailed loss breakdowns or data are not available. 

12. The systematic investigations into secondary loss presented 
in Refs. 18 and 19 are of real interest although the test data is too 
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limited to permit any universal application. An attempt is made to 
define the two-dimensional loss by theory, measure the end loss at 
the blade root and the clearance loss at the tip, and subtract these 
from the total pressure loss measured to arrive at the secondary 
loss. The tests are performed on a particular airfoil shape in 
three solidity arrangements and three flow deflections. The inac¬ 
curacies and errors discussed should help any further investigator 
in this subject. 

13. A comparable breakdown of the losses into components as 
assumed here is offered in Ref. 20. This work by Markov also includes 
the experimental results of some other Russians in the field. With¬ 
out attempting to completely describe the development of equations 
and data presented, the pertinent data for the present consideration 
is extracted. 

The profile loss, or two-dimensional loss which would be caused 
by blades of infinite length is shown on Fig. 10 along with other 
data to be described later. The basis of the Ref. 20 data is a whole 
series of experimental measurements of widely used cascades for 
turbines in Russia. For stators, the profile loss coefficient is 
almost constant at .03 over a wide range of conditions. 

The tip clearance loss is shown on Fig. 9 as represented by two 
plots in Ref. 20. One plot is shown as the loss for a reaction stage 
and one plot is based on experimental data for the efficiency of stages 
with a small degree of reaction. The magnitudes of the degree of 
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reaction are not specified but the trend is correct. Stages with a 
high degree of reaction have a higher loss across the tip, as would 
be expected due to the pressure difference across the rotor from 


entrance to exit. 

The trailing edge loss for cases where the pressure gradient 
is small along the surface of the trailing edges of the blades is 
computed by 



(14) 


Where such a pressure gradient is present the use of boundary layer 
thickness parameters is required. The trailing edge loss coefficient 
according to eq. (14) is shown in Fig. 20. 

The energy lost due to the secondary flow is treated as a 
function of the mass flow rate in the cascade. The efficiency of a 
cascade of blades with finite length is defined as 



(15) 


The loss coefficient is defined as 



(16) 


and so 
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The secondary flow loss coefficient 
used for comparison in this thesis by 


is converted to the 



5^ 5^ 7 (2 ' dim) 

An expression for ^ is developed that is a function of the circu¬ 
lation around the blades, the blade width, and incompressible flow, 
namely 



aa; 


where the coefficient K depends on the expansion of the flow in the 
cascade. The value of K is determined from experimental data to be 


K = 0.0/ 




Thus 


f r = 00/ 


r ”i 

L J 


(17) 


In Fig. 11 the function ^ is shown for the two-dimensional efficiency 
of Ref. 1, 

^7 — a> z = 0.9262 

/(2-dim) 'o 


Values of h/c plotted are chosen arbitrarily but Ref. 20 specifies 
that the derivation holds for blades that are not excessively short. 
There is assumed to be some two-dimensional flow in the channel be¬ 
tween the blades at the center section. Experimental data are shown 
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for short blades, and these data commence at 



■2.0 and go to 


smaller values. So the plot in Fig. 11 is based on eq. (17) and 


is shown for a value of |-j^J ■ 2.0, the lower limit permitted. For 
[ —] ■ 1.33 and 1.25 the mean coefficients are about 107, of the rep¬ 


c \ 


resentative values based on the Vavra equation (3), but these are 
specific data for short blades not arrived at using the equation 
above. The secondary loss as computed by Markov is low compared 
to any other available data. When complete loss coefficients are 
computed in Ref. 20 the values are comparable to any other method 
however, since the tip clearance losses, trailing edge losses, 
blade height corrections, etc., used bring the total loss to the 
proper level. Any use of the Ref. 20 data must take account of 
this fact. 

Of interest also in Ref. 20 is a relation for the discharge 
flow angle of a turbine blade row. The relation is 



(18) 


and it gives more accurate results than the other widely-used 
equation 



(19) 


which does not take account of trailing edge thickness. 
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14. The tip clearance loss as represented by Ref. 21 is also 


shown in Fig. 9. The equation derived is 


* U/-1) + « 



( 20 ) 


whe re 


- .1011 + 4.667 (k/b) 

- -.0422 + 2.790 (k/b) 
E 


h * blade height 


b ■ blade width 


k * clearance 


and the results are shown for two values of (h/b). 

Also given in Ref. 21 is the equation for the minimum induced 
drag due to flow in a clearance gap at the tip of an airfoil (as 
derived by Betz) to be 


2 
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where ^ is a form factor depending on the clearance, blade spacing 
and relative exit angle of the flow. 

Assuming the axial component of velocity to remain constant 
and drawing an arbitrary turbine stage velocity triangle (such that 
all angles are positive), Fig. 12a, the circulation ( P ) can be 
computed by integrating around the closed curve in Fig. 12b to give 



( 22 ) 
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Then the lift per blade is 


L ~~ “"/vie/ ) (23) 


The lift and drag forces exerted on a blade can be seen in Fig. 13 
as well as the axial and peripheral components of the resultant 
force F. It can be seen that for each blade channel 



(24) 



(25) 


The secondary flow loss will be associated with the loss from the 
pressure drop in the axial direction. 

To arrive at a loss coefficient expression for the axial force 
component, it is necessary to define the loss coefficient in terms 
of the pressure loss, as shown in Fig. 14. The loss coefficient is 

defined as 



(26) 


Substituting 
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and for incompressible flow, 


= /^iy^ 


and collecting terms, 



Then from Fig. 13, 






Using the above derived expressions for L and 0 and equating 


(28) 


(29) 




(30) 


gives 

= 7 - i* 

For an arbitrary value of ^(li) ’ ^is secondary loss is plotted in 
Fig. 15. The value of ^ hj can * >e a( *justed to change the general 
inclination of the family of curves shown to any desired value. That 
is, a lower value of K 0 [ , r] raises the curves and vice versa. 
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15. In Ref. 22 an attempt Is made to find a mathematical 
solution to the problem of secondary flows in particular. The 
approximate relations derived are not able to satisfy the boundary 
conditions. One conclusion drawn from the derivations is that the 
secondary flows are a direct consequence of a non-uniform approach 
velocity, and, if the lift is assumed proportional to the square 
of the local velocity, then the secondary flow is independent of 
aspect ratio. These conclusions are counter to the ideas of most 
investigators in the field. 

The tests and measurements made showed odd discrepancies and 
it was felt that rotor-stator induced effects may have been the 
cause. No relation was deduced which is directly useful here. 

16. The ideas and relations presented in Ref. 23 and 24 were 
carried through and refined and elaborated in Ref. 25. All three 
references have the same authors. Therefore, Ref. 25 will be dis¬ 
cussed here as being most useful, and parts of Ref. 23 and 24 will 
will be extracted as necessary. 

By far the most detailed method presently available for deter¬ 
mining loss coefficients is presented in Ref. 25. The loss is 
broken into component parts similar to those assumed at the beglning 
of this discussion and the quantitative data offered are based on 
tests of turbines. 
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The loss coefficients of Ref. 25 are 


Y 



Loss of total head pressure 

Total pressure at blade outlet-static pressure at blade outlet 

(32) 


The loss coefficients used in this thesis are 


f =/ -y 



(33) 


and the necessary conversion is 



For profile losses, two figures are presented, one for nozzle 
blades and one for impulse blades. These are reproduced here as 
Fig. 16a and Fig. 16b. An interpolation method is proposed for rotors 
with reaction and a method of extending the profile loss to conditions 
of other than zero incidence is also shown. 


A number of relations for the loss due to tip clearance are 

presented in Ref. 25 as 



*' Hi) 


These relations are shown in Fig. 9. 


(Stodola, Reaction turbine) (35) 

(Meldahl) (36) 

(507. Reaction turbine) (37) 
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Expressions are also derived for tip clearance loss based on 

and c if the ratio spacing/blade height is small and the gas turning 

JL 

angle is small. For our uses here, the gas turning angles are large. 

The loss relation arrived at is , 

X& - p*/~<.][^ (M) 

where B » 0.5 for blades with clearance and B » 0.25 for blades 
with shrouds. 


Equation (38) is used to find the range of loss coefficients 
for the impulse conditions of cos m 1.00 and ■ 0 (Shown 

in Fig. 9 is the band of values corresponding to o (^ * 0(^ » 60° to 75°. 


The data available to Ainley and Mathieson led them to the con¬ 
clusion that the effect of area ratio is more important to secondary 
flows than is the magnitude of the gas deflection. The area ratio 

used is 


A 2 (annulus area at outlet from blade row) (cos o( 

■■ m .. ■ . ■ . ■■■■■■ 1 ■■■■■■ ..... 

A. (annulus area at inlet to blade row) (cos^^) 

where o( ^ » gas outlet angle measured from the axial and » blade 
inlet angle measured from the axial. The hub ratio (inner diameter/ 
outer diameter) was also felt to exert a strong influence. This led 
to the plot of /Cl + (ID.0D)] versus a parameter ^ used in 

defining the drag coefficient due to the secondary flow as 



( 40 ) 
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This plot is reproduced here as Fig. 17. Ultimately then, a loss 
coefficient for secondary flows is presented as 


K - 






Converting these values of to a plot is presented as Fig. 

for a wide span of inlet and outlet angles. 


( 41 ) 

18 


For trailing edge thickness losses, the assumption is made that 
the previously derived losses are for a blade of trailing edge thick¬ 
ness which is 27. of the blade pitch. For other trailing edge thick¬ 
nesses a correction is applied according to a plot which is reproduced 

here as Fig. 19. 


C- Results of Loss Investigation 

The results of the loss investigation will be discussed according 
to the original breakdown of loss components mentioned at the beginning 
of this section of this thesis. The specific turbine design used in 
the last section of this report will be used as the object of a 
quantitative comparison of the loss relations shown here. 

1. The profile loss estimation method presented by Ainley in 
Ref. 25 is straightforward and complete and should give good data 
for the "conventional" blade shapes from which the data is derived. 

An anomaly presents itself in Fig. 10 however, when a comparison is 
made between the profile loss coefficient data of Ainley and the data 
of Markov. The blade shapes from which the Markov data are derived 
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are very similar to those of Fig. 25. It would appear that such blade 
shapes would have higher profile losses than those based on airfoil 
designs which usually are used for lower flow deflections. The data 
comparison shows this is not the case. 

fa 

Of special interest in the Markov data in Fig. 10 is the fact 
that an optimum solidity of (s/c) ^ .625 to .650 is evident. Higher 
solidities apparently block the flow thereby increasing the profile 
loss, and lower solidities encounter separation effects since the 
blades act more like individual airfoils. The low point in the curve 
corresponds to guiding the flow in a channel, retarding separation 
onset, but not blocking the flow. 

2. The tip clearance losses shown in Fig. 9 seem to give a 
large range of values depending on which curve is chosen. The sit¬ 
uation is perhaps more reasonable when the ratio of radial clearances 
to blade height is restricted to realistic values less than .03. 

The loss estimate from the Ainley equation is obviously outside the 
flow deflection range intended for that relation. The other curves 
for tip clearance loss for reaction blading could be based on dif¬ 
ferent degrees of reaction corresponding to the fan aspect of the 
curves, i.e., high reaction blading would use the highest curve. 

Some reasonable estimate could be made. The only relation for 
impulse blading available is that of Markov. The Meldahl relation 
is the only one which accounts for the trailing vortex inducing a 
loss even though the tip clearance is zero, but the dependence on 
blade aspect ratio appears unreasonably strong for this data as a 
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whole. For (h/b) ■ 1.00, the loss at (k/h) ■ .02 is greater than the 
highest ordinate shown. 

3. The Markov data is all that is available for making an 
estimate of trailing edge thickness loss directly. The plots in 
Fig. 20 indicate the trends which would be expected. 

The Ainley method of modifying the complete loss coefficient 
depending on the relative trailing edge thickness seems an unfor¬ 
tunate way to present the data. The system is workable but indirect. 
It would be more useful for comparison and design purposes to have 
the trailing edge thickness data presented directly. 

4. The situation regarding secondary flow loss coefficients is 
the most confused of all. The general equation by Vavra is not very 
precise but certainly indicates the proper trend if the model pre¬ 
sented in Fig. 7 is correct. None of the other presently available 
methods investigated will even show the proper trends for high flow 
deflections. 

The Ainley relation plotted in Fig. 18 seems to have depended 
too heavily on stator data with an axial entering angle and low flow 
deflections in the rotor. Otherwise, the curves indicate that for 
a given entering angle, increasing the flow deflection decreases the 
loS8. 

In Fig. 15 the relation derived from Meldahl and Betz indicates 
the proper trend for low flow deflections only. Taking the plot for 
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^ • -45° as an example, it can be seen that for flow deflections 
greater than 3 90° the indication is that increasing the flow 

deflection decreases the loss. 

The Markov relation plotted in Fig. 11 suffers the same deficiency 
pointed out for the Ainley and Meldahl curves plus abnormally low 
secondary loss coefficients throughout for reasonable blade aspect 
ratios. The data is of no value unless used specifically with Markov 
relations for all losses. In Ref. 20, when sample loss estimations 
are carried out, the total loss coefficients compare well with the 
measured values. Apparently the profile, tip clearance and trailing 
edge loss coefficients with the blade height corrections used in the 
complete method counter the low secondary loss estimates. 

5. The relations now available can be used to make a quantitative 
comparison of the component losses. The turbine designed in the 
following section will be used as a model to show the comparisons. 

The physical dimensions of the stator and rotor are as listed in 
Table IV and these are used in the relations and plots so far developed. 
The numerical results are shown in Table V for the stator and Table VI 
for the rotor. 
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Table IV 

TURBINE PHYSICAL DIMENSIONS 



Stator 

Rotor 

Entering angle 

0° 

69.60° 

Exiting angle 

75° 

-71.57° 

Blade height, h 

1.165 in 

1.465 in 

Trailing edge thickness, t g 

.036 in 

.025 in 

Blade chord, c 

3.16 in 

1.355 in 

Blade width, b 

1.60 in 

1.375 in 

Spacing, s 

1.80 in 

.9415 in 

Blade thickness, t 

. 69 in 

. 648 5 in 

Tip clearance assumed, (k/h) 

o 

o 

.02 


Table V 

STATOR LOSS COEFFICIENT COMPARISON 


Author 

Profile 

Loss 

Tip 

Clearance 

Loss 

Trailing 

Edge 

Loss 

Secondary 

Loss 

Vavra, Ref.l 




.070 

Ainley, Ref.23 

.067 




Ainley, Ref.25 

.047 



.041 

Markov, Ref.20 

.030 


.032 

<.005 

Meldahl, Ref.21 


.106 


.022 

Stodola, Ref.3 

.040 
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Table VI 

ROTOR LOSS COEFFICIENT COMPARISON 


Author 

Profile 

Loss 

Tip 

Clearance 

Loss 

Trailing 

Edge 

Loss 

Secondary 

Loss 

Vavra, Ref.l 




.231 

Ainley, Ref.23 

~ .12 




Ainley, Ref.25 

.148 

.13 


.168 

Markov, Re f.20 

.068 

.025 

.032 

.025 



* 



Meldahl, Ref.21 


.223 


.2601 

Stodola, Ref.3 


.05* 




* 

reaction 


As can be seen, the component losses presented in Tables V and VI 
can be combined selectively to obtain almost any desired answer. 

The loss investigation is disappointing in that no clear-cut 
answers are provided. The complete Ainley method, or the complete 
Markov method, may produce a reasonable estimate but the correlation 
is poor at best for component losses. 

Any desire to find a way of basing a new design on component 
loss considerations is presently thwarted. A need is evident for 
an orderly testing program to find the magnitude and interaction 
of component losses. 
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IV. Design of Turbine 


As a further investigation into the characteristics of a turbine 
being used for maximum work output, a single stage turbine is designed 
for a slightly supersonic velocity leaving the stator and a relative 
Mach number of 0.8 as seen by the rotor. Then the overall pressure 
ratio is increased and the off-design performance evaluated. The 
number of calculations required are prohibitive for hand calculations 
so the equations necessary for the performance evaluation are also 
programmed in Fortran language for the CDC 1604 Digital Computer. 

A. Loss Coefficient Data 

A set of loss coefficients for the design were provided by 
Vavra from some unpublished data. These loss coefficients are shown 
in Figs. 21, 22 and 23 and will be used for the design and the per¬ 
formance evaluations. 

The abscissa of the plot of stator loss coefficients, M^ g , is 

the Mach number for an isentropic expansion from the total pressure 

at the inlet to the static pressure at the discharge of the blade 

row. This can be seen in the T-s diagram of the expansion through 

a turbine shown in Fig. 24 to be the isentropic velocity V , divided 

tn 

by the acoustic velocity based on T^ ^ g . When reduced, the result is 



( 42 ) 
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To properly determine the flow area through the stator, the 
loss coefficient in Fig. 20 noted "for flow rates" and shown in 
Fig. 24 as *.i is used. This loss coefficient is intended to 
represent the losses encountered by the flow in the channel between 
blades, but is not intended to account for the mixing losses and 
separation losses encountered by the flow leaving the blade row. 
These additional losses are included in the higher loss coefficients 
which are applied to velocity determinations. These loss coef¬ 



ficients are noted "for velocities" in Fig. 21 and as 


Similar considerations are used for the loss coefficients through the 
rotor. * 

The loss coefficients are corrected for blade height according 
to the curves given in Fig. 23. 

B. Flow function 

For steady adiabatic flow the stagnation enthalpy remains con¬ 
stant along a streamline. For stators, the stagnation enthalpy is 



(A3) 


and in Fig. 24, to a different scale, H q would be represented by 
T q and h^ would be represented by T^. A similar situation pertains 
to a relative stagnation enthalpy defined as 



(44) 


for the condition such that the peripheral velocity at the entrance 
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of the rotor equal to the peripheral velocity at the exit of the 
rotor. As shown in Appendix I, the mass flow rate for a given set 
of inlet conditions can be expressed by 


% A > 


V 


aw 

j Hi 1 PJ 


( 45 ) 


where n is the polytropic exponent, P q is the total pressure ahead 
of the stator blade row and p^ is the static pressure behind the 


blade row. 


As shown in Fig. 24, a relative total temperature and relative 
total pressure can be defined as 


7~RI ~ 

T, + 



/ V -£. 


W 


(46) 


(47) 


The mass flow rate through the rotor is then 


sUr' — 





(45a) 


Rearranging the above expressions for flow rate in non-dimensional 
terms defines a flow function. 





( 48 ) 
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and 


p f&_ _ V ZV_ 

Jr-/ 

It will be necessary to investigate whether this flow function can 
be satisfied for each set of conditions imposed across a blade row. 
The choked flow condition corresponds to the maximum value of P/p 
which the blade row can accomodate. The development of the critical 
values of the flow function and pressure ratio can also be found in 
Appendix I to be 






22 

/ 


(49) 


(50) 


A computer program to do these computations is explained in section E 
below. 


C. Geometry and Velocity Triangles 

The final design of the stator and rotor blade rows is as shown 
in Fig. 25 for the mean diameter. The velocity triangles of the 
average flow velocities at each station are as shown in Fig. 26. 

The sign convention chosen is that positive angles and positive 
peripheral velocity components point in the same direction as U, the 
rotor peripheral velocity. 
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It is convenient to define an "exit plane*' and a "discharge 
plane" leaving each blade row as shown in Fig. 25 by the (e) and (d) 
designations. The design of the stator and rotor is such that the 
minimum cross-sectional area of the channel between blades occurs 
at the exit plane (e^ . For maximum work output, the pressure drop 
across the stator is great enough to cause choked flow at the exit 
plane plus an additional expansion of the flow from the exit plane 
to the discharge plane. The conditions for such after-expansions 
are as shown in Appendix II based on the treatment in Ref. 26. 

In Ref. 1 a similar development is shown for a supersonic flow 
entering a blade row. Since the flow leaving the stator is super¬ 
sonic, it is possible under certain conditions for the rotor to see 
a supersonic flow also and proper account must be taken of the shock 
losses involved. 

D. Dimensionless Parameters 

As has already been mentioned, the number of variables to be 
handled is greatly reduced if the mass flow expression is non-dimen- 
sionalized to form a flow function depending only on pressure ratio, 
specific heat ratio and the polytropic exponent of the expansion 
process. 

In Ref. 27 can be found a complete development of dimensionless 
parameters and referred values for minimizing the number of variables 
necessarily handled in presenting the data for a turbine performance 
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analysis. Those used in the present analysis are referred rpm, over¬ 
all pressure ratio, referred flow rate, and a power coefficient. The 
necessary coefficient forms are computed from 


Referred RPM * RPM/ / T 

9 o 

(51) 

Overall pressure ratio ■ P Q /p 

(52) 

Referred flow rate ■ w ~f T /P q 

(53) 

Power coefficient ■ HP/P "/""t 

o / o 

(54) 


The referred flow rate and power coefficient are constant in this 

evaluation. It is useful to refer all the conditions at any point 

in the turbine to inlet conditions. For this reason, velocities are 

carried as V/ V T and W/ and pressures as p/P and P_/P . 

to f o * 0 R 0 

E. Design Computations 

In order to design the stator for choked flow at the exit plane, 
the flow function must be known for the imposed conditions. A com¬ 
puter program was written to compute the flow function for a given 
specific heat ratio, 3* , pressure ratio, P/p, and loss coefficient, . 

As a by-product the program also computes the polytropic loss coef¬ 
ficient , *V This information is presented in Appendix V as Table III, 
Flow Function and Polytropic Loss Coefficients, for W from 1.25 to 
1.40, ^ from 0.0 to 0.25 and P/p from 1.02 to 6.00. The particular 
flow function data to be applied here for ^ * 1.37 are also included 
as Fig. 27. The computer program is included in Appendix III. 
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The following data were assumed given for the design: 


Mass flow rate 


■ 100. lbm/sec 


Total pressure at inlet ■ 1000 psia 
Total temperature at inlet ■ 1460°R 


Specific heat ratio 


- 1.37 


Gas constant, R 


= 421.5 ft lbf/lbm °R 


Mean diameter of stator * 23.0 in 


Mean diameter of rotor 


- 23.0 in 


Design rpra 


- 13000 


The stator blade design was provided by Vavra from unpublished 
data and is as shown in Fig. 25. It is now necessary to select the 
blade height and number of blades so the choked flow condition at 
the exit is satisfied. 

A blade height of 1.16 inches is taken as a first estimate. Then, 
from Fig. 21 the loss coefficient for flow areas at M^ g * 1.00 is 
found to be $ ao m *053. From Fig. 23 the blade height correction 
is 1C 3 0.872 and so 



Using linear interpolation on Fig. 27 for this loss coefficient 


between the maximum values of 



found that 


(f> - .65811 

r cr 



cr 


= 1.86 
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Recalling the flow function equation (48), it is found that 


A 0 

'throat 

A 

THU OAT 



13.82 


13.82 

.65811 


21.0 in 2 


From Fig. 25 the measured minimum distance between blades at 
the exit plane is a * .448 in. and s 3 1.8 in. is given. The number 
of blades is then 




Tfl)^ 

4s 


40.2 


(55) 


In order to prevent harmonics in the wake pattern from the stator 
as seen by the rotor, a prime number of stator blades is chosen. 
In this case, z ■ 41 blades will suffice. For 41 blades, 

A = 7L (23.0) - 1.763 

41 


so all measured values from Fig. 25 must be scaled by the scale 
factor. 

Scale Factor “ • 0.9794 

1 • O 


The blade height must now be 


h . A__21,0_ 

al za 41.0 (.448 x .9794) 


1.165 in 


(56) 


Recomputing 0 CT w ith a new loss coefficient based on this new blade 
height makes no measurable difference in exit area required so the 
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scacor row design is now fixed. To simplify the iteration on a 
computer program was written to find the choked flow values of ^* cr > 
( p /p) cr and 'jfp. This computer program is included in Appendix III 
and the computed data are presented as Table VII in Appendix V. 


The flow conditions at the stator exit plane are determined 
from Fig. 24 as follows: 


Po 


_1_ - .538 

1.86 


= (A) r - .8458 

7 l Po) 


AIL _ (1L _ZL±t)// - £) ..i47o 

T U Z) 

"X • 121.52 

r e I r-i % 

When a check is now made of the Mach number at the exit plane in 
order to check for choked flow it is found that 



- 0.97 


This computed Mach number is the average Mach number of the flow at 
the minimum area. The highest Mach number at this area is equal 
to 1.00 but it is averaged with the Mach number of the flow in the 
boundary layer as accounted for by the loss coefficient. In 
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Appendix IV the derivation of the Mach number for choked flow in a 
polytropic expansion is presented. The "critical average" Mach 
number is found to be 



—/ 


T-/ 


(IV, 9) 


This condition will also be encountered later in the computer program 
for the flow through the rotor and proper account must be taken of 
it there. 

Using the relation from Ref. 20 for flow exit angle and the 


measured blade exit angle from Fig. 25 of ■ 73.9° gives 




an o( f — 


.268 


— (rfji //C&t, °*8/) 

75° 


(18) 



For the rotor, the blade height is permitted to be approximately 
0.30 inches greater than the stator blade height in order to permit 
the streamlines to expand and offer no resistance to flow leaving the 
stator but not have dead areas above and below the useful flow area 
of the rotor blade. The rotor blade height will be taken then as 
h^ - 1.465 in. 

In order to satisfy the requirement for the rotor to see a 
relative Mach number of 0.8, some after-expansion is required. In 
order to evaluate the proper amount of after-expansion, the equations 
in Appendix II were coded for the computer to solve for any applied 
pressure ratio between the exit and discharge planes for a given 'T' , 
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Mach number and blade angle. This computer program Is Included In 

Appendix III and the solution for ■ 1.37 and o (^ • 75° over a 

range of pressure ratios Is presented here as Fig. 28. The computer 

program can be used for any specific heat ratio, blade angle, Mach 

* t ■ 

number and pressure ratio range. 

After some iteration, it is determined that p,/P » 0.460 should 

l o 

give approximately the correct pressure ratio across the expansion 
area from the exit plane to the discharge plane. 


Pe Pe 


.460 

.538 


0.855 


Now from Fig. 28 it is found that &Q( ti 0.2° and ^ .0008. The 

total loss coefficient to be applied is based on a blade height of 
1.165 inches for the stator, ffcer# is ^ ■ .1045 from Fig. 21 and 
■ .838 from Fig. 23 for velocity computations. Therefore 


f, - t(Kt) + £ - • 0M3 


Now the velocity leaving the stator is determined with Fig. 24 
as reference 


T IM , _ 
T 0 

(t -i 

(^o)- 2701 ■ .«• 

477 - 
T 0 

(f-f*) 

(t — ^ w .1892 (.9117) - .1725 

if "J 

m 131.60 

r *W % 
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To find the relative angle and relative velocity of the flow, it is 
necessary to construct the velocity triangle as shown in Fig. 26, 


R 

yr 




7T_ («PM/fr 0 ) J), 

72o. 


W) 


34.18 


75° - 0.2° • 74.8° 


34.5 


- 127.0 


^n/ _ V/ 74.8° 

frTft 

Ym./ „ Yl <Urt, 74.8° 

JT fZ 

RhU = Ru. — ^ . 92.82 

IT fZ fT 

>== » r-T- • 99.07 

/IT /5T 

As a check on the requirement of a relative Mach number of 0.8 as 
seen by the rotor 


M, = WJZ 


99.07 

123.9 


0.80 


To design the rotor, several considerations must be kept in 
mind. The flow deflection should be as great as is practicable for 
maximum work output. The trailing edge thickness will be arbitrarily 
selected as .025 inches, and the spacing and blade design will be 
chosen so the minimum cross section area of the flow passage between 
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blades will occur at the "exit" plane. The rotor blade width is 
selected as 1.375 inches and impulse conditions only will be used 
for unchoked flow, i.e., p^^ ■ p^. 

Recalling that a relative total temperature and relative total 
pressure are defined as 


5/ . JL + . H+^lhLL ■ .9253 



The flow function as applied to the rotor is 


<t> 


f* 

A* • / 


(48a) 


and this expression will be used to determine the exit area of the 
rotor. For the prescribed impulse conditions 



.6955 

.4600 


1.515 


and the value to be used on the abscissa of the plots for rotor 
loss coefficients in Figs. 22 and 23 



34.18 

131.61 


« .2598 
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With this value the rotor loss coefficient for flows is found to be 


^2 • .819 (.192) • .1572 

Now using Fig. 27 for this pressure ratio and loss coefficient and 
using linear interpolation the flow function is 

ft • .58953 

and from eq. (48a) above 

0A 2 • 19.13" 

Aj • 32,44' in^ 

An impulse rotor operating under conditions as these should have 
a ratio 

(*) ao - 7 

and 60 

s • 0.7 (1.375) • 0.962 

For the mean diameter of 23.0 inches, the number of blades on the 
rotor is then 

*2 * " ^ 962 ^ " 75,1 blades 

Since an odd number of blades was chosen for the stator the rotor 
should have an even number of blades to reduce the possibility of 
resonance induced vibrations and the resulting metal fatigue. 


49- 



Arbitrarily choosing 76 blades 



• .9415 in 


and tha distance between blades, a, to obtain the required area of 
32.44 square Inches Is 



- .273 in 


The exiting blade angle can be seen from Fig. 25 to be 



• M - 71570 


The blade design Itself is based on a circular arc contour of 
the concave face. The geometric relations of the design have been 
programmed by Vavra for computer solution. Using this computer 
program and requiring solutions for 74, 76, 78, 80 and 82 blades, 
the printout of the computer solution is presented as Fig. 29. The 
values necessary to draw the blade design are shown in Fig. 30 and 
have been used to draw Fig. 25. 

The solution for 76 blades appears valid and so the other solutions 
are not used. The resulting |^| • 0.691. A check on this spacing can 
be made by Brillngs Rule, an early rule-of-thumb for steam turbines 
which is found in Ref. 3 to be 
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(57) 


If] • 2 - 5 ,ln l fi * vhere /* • 90 • fin 

j^J ■ 2.5 sin 36.86° • 1.5 

. , 6667 


Also Kef. 28 presents a relation for spacing and blade angles 
based on an aerodynamic load coefficient which is found to 

have a nearly constant value over a wide range of turbine designs. 
The relation is 




(58) 


and so, using proper angles as defined in Ref. 28 for this case 

ffl . §4- sln 20 - 40 ° - - .704 

l b l 2-0 8lo 18 , 43 ° sln 38.83° 


The relations by Brlling and Zweifel show at least that the choice 
of 76 blades is reasonable for this design. The relation of Zweifel 
in particular is based on two-dimensional considerations and does not 
account for blade height. Some unpublished test data available to 
Vavra indicates the value of 0.691 is a good one for this size turbine. 

To complete the computations necessary to determine the power 
output of the turbine for a given set of conditions, the velocity 
triangle after the rotor is computed. 
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For no after-expansion, * /^B2 * "71.57 , and finding loss 
coefficients in the manner previously described, 


.460 

.6955 


o \ .2701 „ 


^ 1 + 4 . _ | ** . .926 

T ° ” r 0 \rj 

..o„o 

% It t / 


.8275 


5 


*0 

Wi 

vr 


7e/ 

---f- - .8543 

£ £ 


sihei^. 

I r-t T, 


- 84.39 


Using relations as before with Fig. 26 as reference, it is found 


nr 


45.93 


and the power output, with the proper conversion factors is 


HP 


■ fj (U 1 V ul • U 2 V u2> < 1414 > • 48 • 682 Mp 


(59) 


These computations are verified by the computer program. The first 
page of the printout of the computer program is presented as Fig. 31. 
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v. Computer Program 


A. General Considerations 

The computer program for the evaluation of the turbine performance 
Is Intended to be as general as possible within the framework of 
the turbine design method used. Any single stage turbine so designed 
can be evaluated by entering new physical dimensions, entrance con¬ 
ditions and loss coefficient data. The static pressures to be im¬ 
posed ahead of and behind the rotor must be programmed by the user 
within the main program. The loss data must be of the same format 
as Figs. 21, 22 and 23. 

B. Main Program 

1. The main program is intended to make the major computations 
and comparisons and provide control. The subroutines are provided 
to do the lengthy iterative processes and repetitive calculations. 

The flow charts for the main program and the subroutines are in¬ 
cluded with the program listings found in Appendix III. 

The basic concept of the computer program is as follows: 

a. With the input data given, compute the design condition 
including the after-expansion at stator discharge. 

b. Reduce the static pressure through the impulse rotor in 
increments, computing velocity triangles and power at each pressure 
increment. The shock condition for supersonic flow entering the 
rotor is computed. Also the flow function through the rotor is 
checked for the choked condition. 
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c. When the choked condition through the rotor is attained, 
the static pressure at the rotor discharge is lowered in increments 
and the after-expansion condition computed. The resulting power and 
velocity and thermodynamic conditions are computed for each pressure 
increment. 

d. The process is terminated when the loss coefficient due 
to the rotor after-expansion is reduced to zero. 

The program will also terminate if the loss coefficient due to 
after-expansion out of the stator goes to zero or if the loss coef¬ 
ficient due to an entry shock into the rotor goes to zero. Any con¬ 
tinuation of the program beyond these limiting conditions would lead 
to conditions which violate the Second Law of Thermodynamics, as 
explained in Ref. 1. 

The loss coefficient reaches a maximum and then declines to zero 
or negative values as the imposed pressure ratio across the control 
area is decreased in increments for the after-expansion case. The 
exact reason for this is not known. A check was made to see if the 
axial component of the Mach number of the flow leaving the blade 
row reached M 9 1.00 when the loss coefficient reduced to zero. It 
is found this is not exactly the case but that the axial component 
of the Mach number is less than one. Further work is needed on this 
point. 


2. Some pertinent comments may be useful about the main decision 
points in the main program. The flow chart in Appendix III shows that 
the first check is for the magnitude of the static pressure between 
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the stator and rotor compared to the static pressure in the exit plane 
of the stator. For a stator designed as a converging-diverging nozzle, 
the Mach number at the exit plane may be greater than ■ 1.00. If 
the static pressure between the stator and rotor is greater than the 
static pressure at the stator exit plane, the shock condition is 
computed at the stator exit. If the pressures are equal, there is 
no flow deflection, and if the static pressure between the stator 
and rotor is lower than the exit plane pressure, the after-expansion 
case is computed. 

The next decision is based on the relative Mach number as seen 
by the rotor. For the supersonic case, the shock of the flow in the 
rotor entry is computed. For subsonic flow the component of relative 
velocity in the direction of the rotor blade entering angle is used 
to compute flow conditions. 

The flow function and critical pressure ratios are used to 
determine choked conditions in the rotor. As can be seen in Appendix IV 
the Mach number of the flow at the rotor exit plane is an average 
value, and the critical Mach number varies with the loss coefficient 
as shown there. Since the impulse conditions for the rotor are 
specified, some care is taken by an iteration process to find the 
p^ • P 2 condition for choked flow before the program continues. 

Once the choked condition in the rotor is attained, the con¬ 
ditions forward of the rotor exit plane remain fixed. The Mach number 
of the flow at the rotor exit plane is set at R, ■ 1.00 and the 
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static pressure behind the rotor is lowered in increments. The after- 
expansion behind the rotor is then computed until the limiting con¬ 
dition is reached. 

C. Subroutines 

1. The subroutine called TRNGL computes the velocity triangle 
based on rpm, absolute velocity and direction or based on rpm, 
relative velocity and direction. 

2. The subroutine called AFTER computes the deflection of the 
flow due to shock or expansion of the flow out of a blade row. The 
flow angle is incremented from the blade angle and the pressure ratio 
is computed and compared to the imposed value. The process con¬ 
tinues until the flow deflection is found which corresponds to the 
pressure ratio imposed. The detailed equations are derived in 
Appendix II, based on Ref. 26. 

3. The subroutine for computing the flow function for a 
pressure ratio imposed across the blade row is named CPHI. The 
polytropic loss coefficient corresponding to the design flow 

loss coefficient is determined by iteration. Then the flow function 
to satisfy the imposed conditions, critical flow function, and 
critical pressure ratio are computed from the equations developed 
in Appendix I. The pressure ratio was found to be more sensitive 
in fewer significant figures than the flow function itself so 
decisions in the main program are based on pressure ratio comparisons. 
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4. Subroutine BEFORE is used to compute rotor entry shocks based 
on the relations in Ref. 1 which in turn are based on considerations 
similar to those of Ref. 26 for after-expansion. Since the flow 
deflection is known, an iteration process is not required and the 
computations are straightforward. For cases which are not directly 
solvable by the basic equations, the component of the velocity in 

the direction of the blade angle is taken as the useful resultant. 

5. The CURVE subroutine is written as a means of finding a 
data point on a curve. The loss coefficient curves "for velocities" 
from Figs. 21, 22 and 23 are stored in the computer in the form of 
arrays. The subroutine writes the second-order equation of the 
appropriate three stored data points nearest to the input abscissa, 
and then finds the ordinate on the curve with the input abscissa. 

The subroutine is arranged in such a manner that the stored data 
must have equal abscissa increments and the first abscissa point 
must be zero. This means arbitrarily extending a curve such as 
Fig. 21 to M. 9 0.0 whether the data at low values of M. will be 
used or not. 

D. Results of Computer Program 

1. The turbine design data given in the preceding section of 
this thesis were entered into the computer program and the resulting 
printout for the design condition only is presented as Fig. 31. The 
pressure ratio across the turbine was then increased in increments 
(by lowering p^ and p^ while holding P q constant) and the resulting 
computations are included as Table VIII in Appendix V. The rpm wa s 
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also decreased from 13000 to 9000 in 1000 rpm increments and the 
complete pressure range imposed at each rpm. These results are 
also included in Table VIII. 

2. The non-dimensional parameters previously discussed were 
then plotted in order to make cross plots which ultimately result 
in a portion of the turbine performance map, 

3. The data immediately available from the computer printout 
permits the direct plotting of Fig. 35, Turbine Efficiency vs. 
Pressure Ratio, with the rpm as a parameter. The anticipated gradual 
loss of efficiency with increasing pressure ratio is evident. Of 
interest is the change in slope of the curve at the point where 
rotor after-expansion commences. The indication is that the increase 
in velocity due to after-expansion causes the peripheral component 

of the velocity to increase, even though a deflection of the flow 
toward the axial direction also occurs. Therefore, the rate of loss 
of efficiency is decreased, although the change is very small. A 
reduction in rpm contributes to a large loss in efficiency at any 
given pressure ratio. 

4. In Figs. 36 and 37 the most important indication is as 

follows: Fig. 37 shows that for a fixed design at a pressure ratio 

of (P^/^) * 2.174, the increase in power for an increase in overall 
pressure ratio is direct and marked up to a pressure ratio of about 
4.5 or 5.0. Further large increases in pressure ratio net very 
small gains. In Fig. 36 the corollary is shown for efficiency. The 
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initial increase in pressure ratio and power output causes a certain 
loss of efficiency but the remaining increase in pressure ratio and 
power output causes a drastic loss of efficiency. 

An interesting comparison can now be made with the high head 
coefficient data from the first section of this thesis. In that 
section, each point on each curve represents a different turbine 
design. The conclusion was that if a required high power output 
transcends the need for high efficiency, then a very high head 
coefficient should be used for the turbine design. The initial 
increase of head coefficient from the optimum value causes a large 
loss of efficiency. Once this initial efficiency loss is accepted, 
further increase in head coefficient require little more loss of 
efficiency. Now in this section concerned with the off-design per¬ 
formance of a particular fixed turbine design, the conclusion is: 
additional power can be obtained at a higher pressure ratio than 
design but only in the initial increase of pressure ratio. The 
extremely high loss in efficiency is encountered if the last possible 
increment of power is demanded. 

5. The turbine performance map presented as Fig. 38 serves to 
indicate the pertinent trends for the parameters investigated. The 
lines of constant efficiency form only a portion of the elliptical 
curves they become over the complete range of referred rpm and 
power coefficient. The only conditions investigated in this thesis 
are at high power coefficients. The peak efficiency for this turbine 
would occur in the lower right hand corner of Fig. 38, i.e., in a 
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lower power coefficient range than shown here and in the referred 
rpm range of about 350. 


6. This turbine designed for high head coefficients can now be 
compared to the theoretical values predicted in the first section of 
this thesis. The head coefficient is defined in eq. (1). The 
isentropic velocity C q is given by reference to the of Fig. 24 as 


2^J~C/3 



(60) 


For pressures and temperatures at the design point 

Referring now to Table I at 0{ * 75°, r* =■ 0.0 and k^ s 3 16.0 the 

values determined from the theoretical conditions can be compared 
to the turbine design. The comparison is shown in Table IX in the 
columns Theoretical #1 and Turbine Design. 



The ratio V g/V ^ was arbitrarily selected equal to one for the 

theoretical development. This was not a selected parameter for the 

turbine design and the discrepancy between the two values of V 0 /V , 

m2 ml 

are the basis for the differences seen in Table IX. Since the ratio 

V 0 /V , was not required to be equal to one in the turbine design, 
m2 ml 

the flow is deflected more in the rotor. This results in higher 
absolute values for ^^ and o (Since the flow deflection is 
greater, the peripheral component of the absolute velocity (V is 
greater and the power output and efficiency are increased. Otherwise, 
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the theoretical predictions are good. The V^/V^ value from the 
turbine design can be entered into the equations and the corres¬ 
pondence should be more exact. This was done and the resulting data 
is the column headed Theoretical #2 in Table IX. 


Table IX 

COMPARISON OF THEORETICAL AND TURBINE DESIGN VALUES 



Theoretical #1 

! 

Turbine Design 

r '■ — 

Theoretical #2 

*1 

75° 

74.8° (after-expansion) 

75° 

* 


0 * 

0 

0 

* 

R 2 /R x 

1.0 * 

1 

1.0 

1 

* 

V m2 /V ml 

1.0 * 

.784 

.784 

* 

9 

.95 * 

.955 

.95 

* 


.95 * 

1.00 

.95 

* 

r 

.849 

.852 

.833 


* 

r 

0.0 * 

0.0 

0.0 

* 

k is 

16. * 

16.2 

16.2 

* 

7 

.60111 

.723 

.60358 


2.263 

2.42 

2.031 


V ml /U l 

.984 

.998 

.990 


w 2 /w 1 

.827 

.852 

.812 


4 

69.78° 

-65.30° 

69.60° 

-71.57° 

69.82° 

-70.54° 



-49.17° 

-59.84° 

-57.02° 



selected, not computed 
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7. The curves of Fig. 26 are shown again in Fig. 39 with all 
values referred to the design point. All values are then shown as 

a percent of the design point values. Fig. 39 demonstrates the strong 
effect of rpm on power output and efficiency. At 707o design rpm, it 
is not possible to develop design power output. At about 907* 
design rpm, the design power can be developed at 957o of design 
efficiency. A maximum power of about 1257* of the design value is 
the maximum which can be developed at 907* design rpm. At the design 
rpm, 1257* design power can be developed with a loss of efficiency 
of 107*. When greater power increments are required, the efficiency 
loss becomes increasingly large. A 407* increase in power can be 
developed but at only 72.57* of design efficiency. 

At design rpm, 1257* design power requires an<overall pressure 
ratio of 3.5, and 1407* design power requires an overall pressure 
ratio of 8.5. 

8. As this thesis was being published, it was determined that 
the computer program finds the off-design performance of a turbine 
in which impulse conditions are maintained until choked flow occurs 
in the rotor. A rotor blade height of 1.648 inches is required to 
satisfy such a condition. When this new rotor blade height is 
introduced into the computer program, the only direct influence is 
on the blade height correction factor for loss coefficients. The 
performance computations change is negligibly small. 
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VI. Conclusions and Recommendations 


A. Conclusions 

1. If a requirement for a high shaft power output turbine is 
encountered for which weight, size and space limitations preclude 
having several stages operating at peak efficiency, then a slightly 
higher head coefficient and slightly lower efficiency compromise 
design should not be made. Designing for power outputs higher than 
optimum should be based on extreme power output per stage. The 
initial loss of efficiency is great when a slightly higher than 
optimum head coefficient is considered. Not much more efficiency is 
lost by imposing extreme values of head coefficient. 

2. From the available literature, no straightforward method 

of computing losses is found which correlates well with the available 
data. This appears to be a fruitful area for future research and 
testing. 

a. The available data indicates the tip clearance loss is 
not zero even for impulse conditions and no clearance gap. The 
trailing vortices induced at the ends of the blades introduce 
velocity components which do no useful work. Clearance gaps within 
the annulus boundary layers should give tip clearance loss values 
about the same as zero clearance. As the clearance gap is increased, 
the loss increases. Also the loss increases as the degree of reaction 
is increased. A high degree of reaction means a high pressure drop 
across the blade row and such a pressure differential is conducive 
to high leakage across the blade tips. 
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The loss due to tip clearance is also affected by the relative 
size of the gap to the blade height. The effect of leakage across 
the tips of the blades on the secondary flow will make measurements 
difficult to separate into tip clearance and secondary flow losses. 

b. The profile loss is a function of blade thickness to 
chord ratio, spacing, flow deflection and probably general blade 
shape. This loss offers the most possibilities for simple and 
accurate measurements in a cascade test rig. 

c. The trailing edge thickness losses are probably in¬ 
sensitive below a certain ratio of trailing edge thickness and 
spacing, and then assume a real importance at higher values. These 
losses, too, should be relatively simple to measure. 

d. The secondary loss measurement offers the most difficulty. 
No known theory accurately predicts the loss. Indications are that 
the secondary loss is primarily dependent on the flow turning angle 
and ratio of blade height to the chord or axial width of the blade. 

The tip clearance effect is also of primary importance and will be 
difficult to divorce in the measurements. 

e. Once the four component losses mentioned above have a 
firm basis in measurements, it will be appropriate to consider 
Reynolds number and Mach number effects as they affect the losses in 
turbine blade rows. 

3. The computer program for the turbine performance analysis 
yields results as accurate as the loss data available. The program 
only requires inputs of loss coefficients in the specified form and 
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a design based on the flow function Cp for area determinations. 

The analysis of the off-design performance of this particular 
turbine indicates that large gains in power output can be obtained 
at pressure ratios two or three times the design value if a 157o drop 
in efficiency can be accepted. Attempting to obtain the last 
possible increment of power, however, requires extreme values of 
pressure ratio and another large loss in efficiency. 

4. Further investigation into the after-expansion behind a 
blade row is needed. A cursory examination of the magnitude of the 
axial component of the Mach number of the flow leaving an after¬ 
expansion showed it ot be less than one when the loss coefficient 
decreased to zero. The exact reason for this loss coefficient be¬ 
havior is not known. 

Be Recommendations 

1. That an organized effort be made to measure and define 
component losses in blade rows in a form useful for design purposes. 

a. It would appear that trailing edge thickness effects 
and profile losses could be measured in cascade test rigs with a 
limited number of profiles, solidities and blade heights. Some 
correlation of these measurements could be made with Ainley and 
Markov. 

b. The tip clearance loss should be measured in a rotating 
turbine test stand. The direct effects would not be hard to measure 
if the effect on secondary loss was not present. Since this effect 
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is present, it will require study to obtain quantitative correlation 
and isolate the tip clearance loss from the secondary flow effects. 

c. The secondary loss will probably have to be deduced by 
subtracting the outer component losses from the overall efficiency 
as measured for the complete turbine blade row. No method is pres¬ 
ently known for measuring the secondary loss directly. 

2. Tests on the after-expansion behind a blade row should be 
made. Physical measurements may shed some light on the reason for 
the loss coefficient reaching a peak and then decreasing to negative 
values as the expansion is increased. Further theoretical study is 
also required, especially to consider three-dimensional effects. 

The testing program should include this variable. 
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: FIG. 9> 

TIP CLEARANCE 
LOSS COEFFICIENTS 
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FIG . IZ 

VECTOR DIAGRAM OF FLOW THROUGH AN ARBITRARY 
TURBINE CASCADE ; 





















-79- 


/7<S. /3 

VECTOR. P/ACRAM OF FORCES 
ON ONE BLAOE 





















FIG,. /4- 

H-j. D/AG,RAM FOR THE RELATIVE FLOW 
MA ROWR. 
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Profile-loss coefficients for'conventional section blades at zero 
Incidence, tic = 20 per cent; Re = 2 x 10 5 ; M < 0-6. 
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F-J^/9, Effect of trailing-edge "thickness on blade loss coefficients. 
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F/C. Z3 


EFFECT OF BLADE HEIGHT 
ON LOSS COEFFICIENTS 












FIG. Z4- 

T- s DIAGRAM OF EXPANSION PROCESS 
THROUGH A TURBINE 




































STATION O 


FIG. 25 

STATOR AND ROTOR BLADE DESIGN 
AT THE MEAN DIAMETER 

NUMERICAL VALUES FOUND IN TABLE HZ' 
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VELOCITY TRIANGLES FOR TUFB/NE FLOW 
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FLW FUNCTI 

PRISSURE 

RTIO 


1 J 2 
1 J*» 
1 J6 
1J8 
1.10 
2 

1.14 
1J6 
1J8 
120 
122 
124 
126 
128 
1 JO 
1J2 
1.14 
1 26 

1.8 
1.0 
.2 

1.4 

1.6 

1.8 
1.0 

1.2 

1.4 

1.6 

1.8 
1.0 

1.2 

1.4 

1.6 

1.8 
1.0 

1.2 
1. 4 
1. 6 

1.8 

1.0 

1.2 

1.4 
1.6 

1.3 
1.0 
1.2 
1 .*4 
1.6 
1.9 
2.0 
2.0 
2.0 
2.0 
2.0 
3.0 

3.3 
3.0 
3.0 
3.0 
4.0 

4.3 
4.0 
4.0 
4.0 
5.0 

5.3 
5.4) 
5.0 
5.0 
6.0 









FLOW FUNCTION FOR GAMMA * 1.37 


PRESSURE 

RATIO 

ZETA = 

.000 

.025 

.050 

.075 

1.02 


.19589 

.19340 

.19088 

.18833 

1.04 


.27145 

.26797 

.26444 

.26087 

1.06 


.32588 

.32165 

.31738 

.31305 

1.08 


.36898 

.36415 

.35926 

.35431 

1.10 


.40465 

.39930 

.39389 

.38842 

1.12 


.43496 

.42915 

.42329 

.41735 

1.14 


.46114 

.45493 

.44866 

.44232 

1.16 


.48403 

.47746 

.47081 

.46411 

1.18 


.50423 

.49732 

.49034 

.48330 

1.20 


.52217 

.51496 

.50767 

.50032 

1.22 


.53819 

.53069 

.52312 

.51548 

1.24 


.55256 

.54479 

.53696 

.52906 

1.26 


.56549 

.55748 

.54940 

.54125 

1.28 


.57715 

.56891 

.56060 

.55223 

1.30 


.58769 

.57924 

.57072 

.56213 

1.32 


.59724 

.58858 

.57986 

.57108 

1.34 


.60590 

.59 705 

.58814 

.57917 

1.36 


.61376 

.60473 

.59564 

.58649 

1.38 


.62090 

.61170 

.60244 

.59312 

1.40 


.62739 

•61802 

•60861 

.59913 

1.42 


.63328 

.62377 

.61420 

.60457 

1.44 


•63863 

.62898 

.61926 

.60949 

1.46 


.64349 

.63370 

.62384 

.61394 

1.48 


.64790 

.63797 

.62799 

.61796 

1.50 


.65189 

•64184 

.63174 

.62158 

1.52 


.65551 

.64533 

.63511 

• 62484 

1.54 


.65877 

•64848 

.63815 

.62776 

1.56 


.66171 

.65131 

.64087 

• 63038 

1.58 


.66435 

.65384 

.64330 

.63271 

1.60 


.66671 

•65610 

.64546 

.63478 

1.62 


•66882 

•65812 

.64738 

•63660 

1.64 


.67069 

.65990 

.64907 

.63821 

1.66 


.67234 

.66145 

.65054 

•63960 

1.68 


.67378 

•66281 

.65182 

•64080 

1.70 


.67503 

.66399 

.65291 

•64181 

1.72 


.67611 

•66498 

.65384 

•64266 

1.74 


.67702 

•66582 

.65460 

.64335 

1.76 


.67778 

•66651 

.65522 

.64390 

1.78 


.67839 

.66705 

.65569 

.64432 

1.80 


.67887 

.66746 

.65604 

• 64460 

1.82 


.67922 

.66775 

.65627 

.64477 

1.84 


.67946 

.66793 

.65638 

.64483 

1.86 


.67959 

.66799 

.65639 

.64479 

1.88 


.67961 

.66796 

.65630 

.64464 

1.90 


.67954 

.66783 

.65613 

.64441 

1.92 


.67938 

.66762 

.65586 

• 64410 

1.94 


.67913 

.66732 

.65551 

.64371 

1.96 


.67881 

.66695 

.65509 

.64324 

1.98 


.67841 

.66650 

.65460 

.64271 

2.00 


.67794 

.66598 

.65404 

•64211 

2.20 


.67027 

.65794 

.64564 

.63338 

2.40 


.65909 

.64649 

.63395 

•62146 

2.60 


•64613 

.63334 

.62063 

•60800 

2.80 


.63237 

.61945 

.60663 

.59391 

3.00 


.61838 

.60538 

.59249 

.57972 

3.20 


•60450 

.59144 

.57852 

.56573 

3.40 


.59092 

.57783 

.56489 

.55211 

3.60 


.57775 

.56465 

.55172 

.53896 

3.80 


.56505 

.55195 

.53904 

.52631 

4.00 


.55284 

.53976 

.52688 

.51419 

4.20 


.54113 

.52808 

.51524 

.50260 

4.40 


.52992 

.51690 

.50410 

.49152 

4.60 


.51918 

•50620 

.49346 

• 48094 

4.80 


.50891 

.49597 

•48328 

.47083 

5.00 


.49908 

•48618 

.47355 

.46117 

5.20 


.48966 

.47682 

.46425 

.45193 

5.40 


•48065 

.46786 

.45534 

•44310 

5.60 


.47200 

.45927 

.44682 

.43465 

5.80 


.46371 

.45103 

.43865 

.42655 

6.00 


.45576 

.44313 

•43082 

.41879 


F/G. 27 


.100 

.125 

.150 

.175 

.18574 

.18312 

.18046 

. 17776 

.2572 5 

.25358 

.24987 

•24610 

.30867 

.30423 

.29973 

.29518 

.34931 

•34424 

.33911 

.33392 

.38299 

.37729 

.37162 

.36588 

.41136 

.40529 

.39915 

.39293 

.43591 

.42942 

.42286 

.41623 

.45732 

.45047 

.44354 

.43652 

.47618 

.46899 

.46171 

.45436 

.49289 

.48539 

.47780 

.47013 

.50777 

.49999 

.49212 

.48416 

.52109 

.51304 

.50491 

.49669 

.53304 

.52474 

.51637 

.50791 

.54378 

.53526 

•52666 

.51797 

.55347 

.54474 

.53592 

.52703 

.56222 

.55329 

.54427 

.53518 

.57012 

.56100 

.55181 

.54253 

.57727 

.56798 

.55861 

.54916 

.58373 

.57428 

.56475 

.55513 

.58959 

.57997 

.57029 

.56052 

.59437 

.58512 

.57528 

.56 538 

.5996 6 

.58976 

.57979 

.56975 

.60397 

.59394 

.58385 

.57368 

.60736 

.59771 

.58750 

.57721 

.61137 

•60110 

.59076 

•58036 

.61451 

.60413 

.59369 

.58318 

.61733 

•60684 

.59630 

.58568 

.61934 

.60926 

.59861 

.58790 

•62208 

.61139 

.60065 

.58985 

.62405 

•61328 

.60245 

.59156 

.62579 

.61492 

•60401 

.59304 

.62730 

.61636 

.60537 

.59432 

•62862 

.61759 

.60652 

.59540 

.62973 

.61864 

.60749 

.59630 

•63068 

.61951 

.60829 

.59703 

.63146 

•62022 

.60894 

.59761 

.63208 

.62078 

.60943 

.59805 

.63257 

.62120 

.60980 

.59835 

.63292 

.62149 

.61003 

.59853 

.63314 

•62166 

.61015 

.59860 

.63326 

.62172 

.61015 

.59855 

.63326 

.62167 

.61005 

.59841 

.63316 

•62152 

.60986 

.59817 

.63297 

.62129 

.60958 

.59784 

.63269 

•62096 

.60921 

.59744 

.63234 

.62056 

.60877 

.59695 

.63190 

.62008 

.60825 

.59640 

.63139 

.61953 

.60766 

.59578 

.63031 

.61892 

.60701 

.59510 

• 63018 

.61825 

.60631 

.59436 

.62114 

.60892 

.59672 

.58453 

.60902 

.59663 

.58428 

.57195 

•59544 

.58295 

.57050 

.55812 

•58128 

.56874 

.55626 

.54387 

.56706 

.55450 

.54204 

.52966 

.55307 

.54052 

.52809 

.51576 

.53947 

•52696 

.51457 

• 50230 

.52635 

.51389 

.50157 

.48937 

.51375 

.50135 

.48910 

•47700 

.50169 

.48936 

.47719 

•46518 

. 49016 

.47791 

•46582 

• 45390 

.47915 

.46698 

.45498 

.44317 

.46864 

.45655 

•44465 

.43293 

•45861 

•44660 

.43479 

.42318 

.44902 

.43710 

.42539 

•41388 

.43987 

.42803 

.41642 

• 40501 

•43111 

.41937 

.40784 

• 39654 

.42274 

•41108 

.39965 

•38845 

.41472 

•40315 

.39182 

•§8071 

.40704 

.39555 

.38431 

.37331 
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.200 


.225 


.250 


.17503 

.24228 

.29055 

.32865 

.36006 

.38663 

.40950 

.42942 

.44691 

.46238 

.47612 

.48838 

.49935 

.50920 

.51804 

.52600 

.53316 

.53962 

.54544 

.55068 

.55539 

.55963 

.56343 

.56684 

•56988 

.57259 

.57500 

.57712 

.57899 

.58061 

.58201 

.58321 

.58422 

.58505 

.58572 

.58624 

.58662 

.58687 

.58699 

.58700 

.58691 

.58673 

•58644 

•58608 

.58563 

.58511 

.58452 

.58387 

.58316 

.58239 

.57235 

.55965 

.54578 

.53153 

.51736 

.50351 

.49014 

.47730 

.46502 

.45331 

.44214 

.43151 

.42139 

.41175 

•40256 

.39379 



.17225 

.23840 

.28587 

.32330 

.35416 

.38025 

.40270 

.42223 

.43938 

.45453 

.46799 

.47998 

.49071 

.50033 

.50896 

.51673 

.52371 

.52999 

.53565 

• 54074 

.54532 

.54942 

.55310 

.55639 

.55933 

.56193 

.56424 

.56627 

.56805 

.56959 

.57091 

.57204 

.57297 

.57374 

.57435 

.57481 

.57513 

.57532 

.57540 

.57537 

.57523 

.57500 

.57468 

.57427 

.57379 

.57324 

.57262 

.57193 

.57119 

.57040 

•56016 

.54737 

.53347 

.51925 

.50513 

.49136 

.47808 

.46534 

.45317 

.44157 

.43052 

.42001 

.41001 

•40049 

.39141 



. 16942 
.23446 
.28111 
.31788 
.34817 
.37378 
.39579 
.41494 
.43174 
.44658 
.45975 
.47148 
.48196 
.49136 
.49978 
.50735 
.51416 
.52027 
.52577 
.53071 
.53515 
.53912 
.54268 
.54586 
.54868 
.55119 
.55340 
.55534 
.55703 
.55849 
.55974 
.56079 
.56166 
.56237 
.56291 
.56332 
.56359 
.56373 
.56376 
.56368 
.56350 
.56323 
.56287 
.56242 
.56191 
.56133 
.56067 
.55996 
.55919 
.55838 
.54796 
.53511 
.52120 
.50701 
.49297 
.47929 
.46611 
.45349 
•44144 
.42996 
.41904 
.40866 
♦39878 
.38938 
.38044 
17191 

m 

14863 
14156 




MACH 

ALFA 

GAMMA 



1.000 

75 

T.370 

• 


PRAT 

PRAT 1 

VELRAT1 

CALFA 

ZA 

1.0000 

1.0000 

1.0C00 

-.0000 

.0000 

.9800 

.9800 

1.0146 

.0039 

• 0000 

.9600 

.9600 

1.0292 

.0157 

.0000 

• 9400 

.9401 

1.0438 

.0355 

.0001 

.9200 

.9201 

1.0583 

.0636 

.0001 

.9000 

.9001 

1.0729 

.1002 

.0003 

.8800 

.8799 

1.0876 

• 1460 

.0005 

• 8600 

.8601 

1.1C21 

.2001 

.0008 

• 8400 

• 8401 

1.1167 

.2641 

• 0011 

.8200 

.8201 

1.1313 

.3380 

• 0016 

.8000 

.7999 

1.1461 

.4230 

.0021 

.7800 

.7800 

1.1606 

.5180 

.0028 

.7600 

.7599 

1.1753 

.6250 

.0036 

.7400 

.7399 

1.1900 

.7440 

.0045 

.7200 

.7200 

1.2045 

.8750 

.0055 

.7000 

.7001 

1.2191 

1.0200 . 

.0067 

• 6800 

.6800 

1.2338 

1.1810 

.0080 

• 6600 

.6600 

1.2485 

1.3570 

.0095 

.6400 

• 6400 

1.2632 

1.5500 

.0112 

.6200 

.6199 

1.2780 

1.7620 

.0130. 

• 6000 

• 6000 

1.2928 

1.9930 

.0150 

.5800 

.5799 

1.3076 

2.2460 

.0171 

.5600 

.5599 

1.3225 

2.5230 

.0194 

.5400 

.5399 

1.3374 

2.8250 

.0219 

.5200 

.5199 

1.3525 

3.1560 

.0246 

.5000 

.4999 

1.3676 

3.5180 

.0274 

TIME, 1 

MINUTES 

ANO 36 SECONOS 



FIG. 23 

AFTER-EXPANS/ON AT K-/37AND & t =75° 
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BLADING FOR ROTOR DESIGN 


DIAMETER 

9 

23.0000 

23.0000 

23.0000 

NUMBER CF BLADES 

> * 

. 74 

76 

78 

INLET ANGLE 

S 

20.4000 

20.4000 

20.4000 

EXIT ANGLE 

S 

18.4300 

18.4300 

18.4300 

AXIAL WIDTH 

• S 

1.3750 

1.3750, 

1 .3750 

SPACING 

S- 

.9764 

.9507 

.9264 

EDGE THICKNESS 

S 

.0250, 

.0250 

.0250 

R 

S 

.7291 

.7291 

.7291 

XI 

e 

.254 1 

.2541 

. .2541 

X2 

k 

.2305 

.2305 

.2305 

Y 1 

k 

.6833 

.6833 

.6833 

Y2 

s 

.6917 

.6917 

.6917 

X 3 

s 

.3259 

. 3259 

.3259 

X4 

s 

.3096 

.3096 

.3096 

R$ 

s 

.4276 

.4362 

.4443 

X5 

s 

.9821 

.9562 

.9317 

Y3 

s 

.7002 

.6998 

.6994 

YU 

s 

.6748 

• .6752 

.6756 

X 6 

s 

2.2681 

2.2681 . 

2.2681 

Y5 

s 

.7223 

.7223 

.7223 

Y6 

s 

.6527 

.6527 

,6527 

AREA 

s 

.661387 

.649473 

.637951 

DELTA 

s 

.9850 

.9850 

.9850 

CGLXS) 

s 

.9434 

.9361 

.9292 

CGIX7) 

s 

.9433 

.9360 

.9290 

CGIY7) 

s 

.6996 

.6994 

.6993 

CGIY8) 

s 

.6754 

.6756 

.6757 

I(X-X) 

- • 

06302621 . 

06249739 . 

06196594 

IIY-Y) 

s • 

03595965 . 

03399491 . 

03219183 

X I MAX ) 1 

s 

.6378. 

.6305 

.6235 

Y(MAX) 1 

s 

.6619 

.6619 

.6619 

X I MAX ) 2 

s 

.6259 

.•6186 

.61 16 

Y I MAX) 2 

3 

.6876 

.6876 

.6876 

X(MAX)3 

3 

.4664 

.4564 

.4470 


AlHMA/J * 

TIME, 0 MINUTES AND 33 SECONDS 


23.CC0C 

8C 

20.40C0 
18.4300 
1.375C 
.9032 
.0250 
.7291 
.25*41 
.2305 
.6833 
.6917 
.3259 
-.3096 
.*4521 
.9084 
.6990 
.6760 
2.2681 
.7223 
.6527 
•6266CQ 
.9850 
.9225 
.9223 
.6992 
.6758 
06143323 
03053410 
.6168 
.6619 
.6049 
.6876 
.4381 


2 3.0000 
82 

20.4000 
18.4300 
1.3750 
.8812 
.0250 
.7291 
.2541 
.2305 
.6833 
.6917 
.3259 
.3096 
.4594 
• .8863 

.6986 
.6764 
2.2681 
.7223 
.6527 
.616030 
.9850 
.9160 
.9159 
.6991 
.6759 
C6090044 
C290C724 
.6104 
.6619 
.5985 
.6076 
.4298 


Ft G. Z9 


! 



F/G. 30 

DIMENSIONS OF ROTOR BLADES 
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TURBINE PERFORMANCE 


INPUT VALUES AND .CONSTANTS COHPUTEO 

FLOW RATE * 100.0 LBM/SEC PO » 1000. PSIA 

“ 11 GAMMA » 1.f“ 


REFERRED FLOW RATE 

Cl * 3.6195 

316.9026 
.000010 
136.3049 


3.1 


GAMMA 


,37 


TO * 
R * 


1460. 3 
421.50 


C4 


EXP1 « .2701 
EXP2 * 3.7027 
EXP3 » 2.7027 
RTO = 38.21 


RPM 


13000. 


PHYSICAL PARAMETERS 

REFERRED RPM 


MEAN DIA AT ROTOR INLET 
MEAN DIA AT ROTOR EXIT 


23.00 IN 
_ 23. 00 IN 

STATOR DISCHARGE ANGLE 
ROTOR INLET ANGLE 
ROTOR DISCHARGE ANGLE 


STATOR BLADE HEIGHT 
ROTOR BLADE HEIGHT 
> 75.00 DEG 

‘ 69.60 DEG 

• -71.57 DEG 


LOSS COEFFICIENT THROUGH ROTOR (FLOWS) 
STATOR EXIT PLANE 


.1572 


340.226 

* 1.165 IN 
> 1.465 IN 


3B 

• < 
om 

VE/RTO 

121.515 

PE/PO 

.5380 

TE/ TO 
.85 30 


. 

STATOR 

DISCHARGE 

PLANE 




MIVll 

1.06 

V1/RT0 

131.612 

P1/P0 

.4600 

T 1/TO 
.82 75 

ALFA1 

74.78 

ZETA1 

.0883 

M(W1) 
.80 

W1/RT0 

99.071 

PR 1 / PO 
.6955 

TRI/TO 
.92 53 

BETA1 

69.60 


IN ROTOR ENTRY 





M(W1) 

• 80 

W1/RT0 

99.071 

PR1/P0 

.6955 

TRI/TO 

.9253 

BETA1 

69.60 


FLOW FUNCTION COMPUTATIONS 

PHI COMPUTED * .58953 
PRESSURE RATIO IMPOSED » 1.512 
ROTOR EXIT PLANE 

CRITICAL 

PHI CRITICAL 
PRESSURE RATIO 

Mt WE) 
.73 

WE/RTO 

90.952 

PE2/P0 
• 4600 

TE2/T0 

.8*29 



ROTOR 

DISCHARGE 

PLANE 




Ml W2) 
.67 

W2/RT0 

84.389 

P2/P0 

.4600 

T2/T0 

.8543 

BETA2 

-71.57 

ZETA2 

.2744 

M(V2) 

.42 

V2/RT0 

53.105 

P2/P0 

.4600 

T2/T0 
• 8S43 

ALFA2 

-59.84 



•60837 

1.827 


RESULTANT POWER AND EFFICIENCY 

VU2CFT/SEC) P2(P , SIAI 
-1754.49 460. 

OVERALL PRESSURE RATIO > 2.174 HP « 48682. 


VUItFT/SEC) 
4852.61 


P2TJPSIA) 

519* 


T2 

1247. 


T2T 

1288. 


EFFICIENCY » .723 


FK*. 3! 



FIG. 3Z 

A'a DIAGRAM FOR POLYmoP/C EXPANSION 
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FIG. 34 

J_ DIAGRAM PDF AFTER-FXPANS/ON PROCESS 
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Appendix I. Development of Flow Function Relations 


In Fig. 32 is shown the adiabatic expansion process in a 

turbine from state 1 to state 2. 

The fluid proceeds from state 1 to state 2 1 and the isentropic 

enthalpy drop is represented by process (1-2!^). The fluid proceeds 

from state 2* to state 2 and the isentropic enthalpy drop is process 

(2*-2 ). For the complete process the isentropic enthalpy drop is 

represented by process 1-2^. However, the sura of processes (1-21^) 

* 

+ (2*-2 ) is greater than process (1-2^) due to the so-called 
M reheat factor 11 . The frictional effects raise the temperature, and 
therefore the enthalpy, of the fluid as it proceeds from state 1 to 
state 2*. This increase of enthalpy is then available to do work. 

As a result, it is necessary to make a distinction between the 
loss coefficient ^ for the overall process of Fig. 14 and the loss 
coefficient f p for the polytropic process from pressure (p) to 
(p-dp) on Fig. 32. 

The polytropic loss coefficient is 

■ 1 - ? P d.i) 


Looking at the expansion process in more detail in Fig. 32, 
the efficiency can also be seen to be 





( 1 . 2 ) 
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From this the relation along the polytropic line is 




% 

If-/ 

/ -{+-4) r 

= X 

3W_ 


( -f J 


V p 


and 






or 


*j-r = % ^ 


r <■ v t 

Since, for an adiabatic process 


'TV ___ v 

a constant, then / oC y* 

and for an isentropic process 

't '^r 

■ constant, then / cC. 

Comparing these exponents, there is 
sr\ — / 




/TV 


- 7 ,(¥) - 


and 


/TL — 


r 


~f p (v-0 ■+ / 


(1,3) 


(1,4) 


(1,5) 


( 1 , 6 ) 


(1,7) 
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To find the relation between the overall loss coefficient and the 
polytropic loss coefficient, eq. (1,5) can be written between two 
points on the polytropic expansion line as 


51 = 

T 


5 

r 


-3 V-! 
'F V 


and 


%-T. r B 



\ Xzf 


v-T 

i-(ri 

7 V 

= 1 T ° 

'-PX* 

U/ 


/ 

UJ 


( 1 , 8 ) 


Solving equation (1,8) for and using 

the reduced result is 

f- 


-r-i 
7 * \ ~ 




K 


' - (?)* 

To define a flow function for the mass flow through the blade 

row, the expression for mass flow rate must be expressed in terms of 
the pressure ratio and the polytropic exponent. The mass flow is 


(1,9) 


•— a 2 v z 


sUS-* = 


"z 


( 1 , 10 ) 


-106- 












and following the procedure outlined in Appendix IV for finding the 
velocity corresponding to a pressure drop the result is 




2 


zr $Z~r~ 0 

r-/ 



(iv,4) 


Along the polytropic line, for an adiabatic process 

/ru 



* constant 


or 






/ 





(i,u) 


When equations (IV,4) and (1,11) are substituted into (1,10) the 
result can be arranged in non-dimensional form to define a flow 
function (j) to be 



JT 



( 1 , 12 ) 


This flow function reaches a maximum at the choked condition for the 
blade row. To find the corresponding pressure ratio for choked flow, 
the derivative of the flow function with respect to the pressure ratio 
will be set equal to zero and the result solved for the pressure ratio. 
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Then 




(tf) 


^Oi 




~/L 

/TV 


/ 


+ 


/n + / 


^>2_±i 

/n 


-/ 


A. 


which reduces to 


Z* 


/y^r 




*w 




/TL sn 


/-** 


When set equal to zero, the result is 


/l 


sn — / 
sru 


= -(—) 
srt Uf// 


and 


P ) 


z 

/n+t 




(1,13) 


When the relation for the critical pressure ratio, eq. (1,13), 
is substituted into eq. (1,12) the critical flow function is obtained 
as 



(1,14) 
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Appendix II. After-expansion from a Blade Row 


To find relations to express the expansion of a sonic or super¬ 
sonic flow out of a blade row, the following model is employed. The 
cascade is assumed to have infinitely many blades of zero thickness 
arranged as in Fig. 33 below. 



Fig. 33 Relations for After-expansion from a Cascade 

A unit height perpendicular to the page is assumed. Also the uniform 
adiabatic flow is assumed to be frictionless at so that R is 

perpendicular to . The control a’rea under consideration is 

bounded by the surfaces (jd) , (e) and (d) , i.e. , the "wall", the 
"exit" plane and the "discharge" plane. V g is sonic or supersonic. 
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The continuity equation is used to find a relation between the 


densities, 


,4 ^ ^ =fit £ C^.(k-Ac<)^ = srn s 

— <d- 'C&Z o( 


and so 


_ (Vt \ C^io( 

/* 


Yt) £#xl(c/-A°<) 


(II.1) 


The momentum equation, Ref. 1, ignoring shearing stresses is 

^ = ^-sYc -^s Yt -A&- A 7^ A ' 


(II.2) 


where 


and so 


k = 


& - 


Jtv. 

A 


- x l^ yC&H A o( — g Ao< 


/n, — 


A 

— A, 


/** — 


A 

— sO 


A 

sCAH C< — J oC 


£ -«n s \/^(a Ao( — J £4* &Aj 4‘U'^ “6- ~^£ ttt o(-j ) 


Of particular interest are the relations in the i direction 

O - /f>1^ 1^ ^ <U^ 4o( -h jbg ~6~ ~ x2- o( 
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and using 


^5 ~ /'jt 'JL 


7* 


k: 


reduces to 

a ^ 

£ L 


Yl-Y* /ten. A c< 


-h / - 

A 


= o 


Since 


= ±- fi = 

f ir 


j/ 

/CL 




/ - V/ /Ott- At/ 

(4 


& 

A 


/ 


(II.3) 


Since the process is assumed adiabatic, the total enthalpy remains 
constant through the control area, 

-i x +£ = 4 + yL 
£ ^ 

and 

-/ = CpT = = J_ JL 

J~ x-i ' r-/ p 

and substituting 

as = n.ik.+ y£ 

*-> A r-/ fij. 

' (11,4) 
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Rearranging and substituting in eq. (11,1) and using 



results in 



/(UXL o( 


/ + 



(11,5) 


Then equating eqs. (11,3) and (11,5) and collecting terms results 
in a quadratic equation in velocity ratio, 



The loss coefficient associated with the flow through the 
control area can be defined to be 






based on the relations in Fig. 34. Using perfect gas relations and 
Fig. 34 as a reference, the following relations are obtained for the 
loss coefficient. 
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The Mach number of flow behind after-expansion is then, using Fig. 34 
as reference for relations, 


^ = (M d-Y—) 

“x 1 vj[ ^){*u) 
H - = J£l 

^L_L— - ^L-i- = ^ ' 

v-i af r-y ^ ^ 



Further 


/ -h rd. M z 








a. 


^ -r-/ 






— 








/ V- 4W /** 
£ -*■ 



( 11 , 9 ) 
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Substitute into (11,9) above 



This computation is included in the computer program but is not 
printed out. It is inserted for ready use if a check is needed. 


( 11 , 10 ) 
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Appendix III. 


MAIN PROGRAM FOR TURBINE PERFORMANCE ANALYSIS 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

SUBROUTINE AFTER 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

SUBROUTINE TRNGL 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

SUBROUTINE BEFORE 
FLOW CHART 
PROGRAM NAMES 
PROGRAM LISTING 

SUBROUTINE CURVE 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

SUBROUTINE CPHI 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

PROGRAMS TURPLOT AND TURBINE (HIGH HEAD COEFFICIENT DATA) 
FORTRAN NAMES 
PROGRAM TURPLOT LISTING 
PROGRAM TURBINE LISTING 

PROGRAM FANNO 3 (FLOW FUNCTION AND ZETA POLYTROPIC) 
FORTRAN NAMES 
PROGRAM LISTING 

PROGRAM FANNO 4 (CHOKED FLOW VALUES OF FANNO 3) 

FORTRAN NAMES 
PROGRAM LISTING 

PROGRAM FOR AFTER-EXPANSION COMPUTATIONS 
FORTRAN NAMES SAME AS "AFTER” ABOVE 
PROGRAM LISTING 
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PROGRAM PERFORM 



FORTRAN NAMES, 

EQUIVALENT SYMBOLS AND DEFINITIONS 

AE2 

A e2 

Area at exit plane of rotor 

AJ 

J 

Mechanical equivalent of heat 

ALFABl 

*.i 

Undeflected flow angle of absolute velocity 
leaving stator 

ALFAl 


Flow angle of absolute velocity leaving 
stator 

ALFA2 

o ( 2 

Flow angle of absolute velocity leaving 
rotor 

AME 

M 

e 

Mach number of flow at exit plane of stator 

AMIS 

M. 

is 

Mach number for isentropic expansion 

AMV1 

*vi 

Mach number of absolute flow leaving stator 

AMV2 

^2 

Mach number of absolute flow leaving rotor 

AMW1 

“wi 

Relative Mach number of flow entering rotor 

AMW2 

^2 

Relative Mach number of flow leaving rotor 

AMWE 

^We 

Relative Mach number of flow at exit plane 
of rotor 

AMWEX 


Average Mach number of flow at exit plane 
of rotor for choked flow 

BETAl 

h 

Flow angle of relative flow entering rotor 

BETA2 

BETAB1 

fh 

y^Bl 

Flow angle of relative flow leaving rotor 

Relative flow angle for zero incidence 
entering rotor 

BETAB2 

y^B2 

Undeflected relative flow angle leaving 
rotor 

CP 

C 

n 

Specific heat at constant pressure 
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Cl 


(Constant #1) = V R/g 

C2 


(Constant #2) - V(2 IT gR) / (tf -1) 

C3 


(Constant #3) = ( -1)/(2 T gR) 

C4 


(Constant #4) = V^gR 

DM1 

°ml 

Mean diameter of blade at rotor entrance 

DM2 

° m2 

Mean diameter of blade at rotor exit 

DR PM 

^RPM 

Increment of rpm 

DTE 

(AW 

Temperature drop from inlet to exit plane 
of stator 

DTF 

Al/T 

O 

Final temperature difference from total to 
static temperature after rotor 

DTR 

<T R1- T 1> 

Temperature drop between relative total 
temperature ahead of rotor and static 
temperature ahead of rotor 

DTRAT 

(iT/T ) 
o 

Temperature drop from inlet to discharge 
plane of stator 

DTE 2 

(ATW 

Temperature drop from inlet to exit plane 
of rotor 

DT2 

(At 2 /t o ) 

Temperature drop from inlet to discharge 
plane of rotor 

EFF 

7 

Efficiency based on total pressure after 
ro tor 

EXP1 


Exponent #1 = (#-1)/^ 

EXP2 


Exponent #2 = V/(Y-1) 

EXP3 


Exponent #3 = l/( T~l) 

G 

s 

Acceleration of gravity 

GAM 


Specific heat ratio 

HKR 

*SiR 

Blade height correction for rotor loss 
coefficient 

HKS 


Blade height correction for stator loss 
coefficient 
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Horsepower 


HP 


HR2 

^2 

Rotor blade height 

HSl 

h si 

Stator blade height 

NMAX 


Number of data points stored for a 
particular loss coefficient curve 

PE 

< p e /P o> 

Static pressure at exit plane of stator 

PE 2 

< p e2 /P o> 

Static pressure at exit plane of rotor 

PHIC 

A 

Flow function computed with input If and ^ 

PHICR 

At 

Critical flow function for choked flow in 
rotor 

PPCR 

( p / p >cr 

Critical pressure ratio for choked flow in 
rotor 

PO 

P 

o 

Total pressure ahead of stator 

PPRAT 

P Rl /p 2 

Pressure ratio through rotor 

PRAT 

p l /p e 

Pressure ratio between stator exit plane 
and stator discharge plane 

PRATO 

P c /p 2 

Overall pressure ratio 

PRP 

P Rl /p l 

Relative pressure ratio 

PR1 

P R1 /P 0 

Relative total pressure ahead of rotor 

PRIY 

P R1* /P 0 

Relative total pressure behind shock in 
rotor entry 

PI 

P l /P o 

Static pressure between stator and rotor 

P2 

p 2 /P o 

Static pressure after rotor 

P2T 

p 2 

Static pressure after rotor 

R 


Gas constant 

RFR 

W V T /P 
/ o o 

Referred flow rate 

RPM 

rpm/ y T q 

Referred rpm 
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RPMA 


Actual value of rpm 

RTO 

f7" 

' o 


TE 

T /T 
e o 

Temperature at exit plane of stator 

TEIS 

T /T 

Isentropic temperature at exit plane of 


e is ° 

stator 

TE2 

T 0 /T 
e2 o 

Temperature at exit plane of rotor 

TI 


True interval between data points on loss 
coefficient curves 

TO 

T 

o 

Total temperature ahead of stator 

Tl 

T./T 

1 o 

Static temperature between stator and rotor 

T2 

T 2 /T o 

Static temperature after rotor 

T2T 

T 2 

Static temperature after rotor 

TP2 

P 2 

Total pressure after rotor 

TRAT 

(T ll s /T o> 

Isentropic temperature ratio through stator 

TRATR 

(T 2 1s /T 81 ) 

Isentropic temperature ratio through rotor 

TR1 

T R1 /T o 

Relative total temperature ahead of rotor 

TT2 

T T2 

Total temperature after rotor 

T1Y 

T * 

1 

Static temperature behind shock in rotor 
entry 

U1 

v/x 

Peripheral speed of blade at mean diameter 
at rotor entrance 

U2 

U 2 / V T „ 

Peripheral speed of blade at mean diameter 
at rotor discharge 

UV 

i—4 

> 

r—4 

Abscissa for rotor loss coefficient data 

VE 

V / V T 
e / o 

Absolute velocity of flow at exit plane of 
stator 

VM 

V / V~T 

m ' o 

Meridional component of absolute velocity 

vul 

V . / T 

ul / o 

Peripheral component of absolute velocity 
leaving stator 
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VU2 

V ./ /T" 

u2 ' o 

Peripheral component of absolute velocity 
leaving rotor 

VI 

v i 7 

Absolute flow velocity leaving stator 

V2 

v-/ y t 

2 » o 

Absolute flow velocity leaving rotor 

WDOT 

w 

Mass flow rate 

WE 

V 

Relative velocity at exit plane of rotor 

WM 

w / /T" 

m ' o 

Meridional component of relative velocity 

QUl 


Peripheral component of relative velocity 
entering rotor 

WU2 

^ -n: 

Peripheral component of relative velocity 
leaving rotor 

W1 

v f*r 

Relative flow velocity ahead of rotor 

WlY 

w */ y t 

1 * o 

Relative flow velocity behind entry shock 
in rotor 

W2 

W 2 / AT 

Relative flow velocity behind rotor 

YY 


Ordinate of loss coefficient curves 

ZET 


Loss coefficient for velocity through rotor 
including after-expansion 

ZETAA 

/ A 

After-expansion loss coefficient 

ZETAB 


Loss coefficient due to rotor entry shock 

ZETAS 

fs 

Stator loss coefficient 

ZETAA1 

fa! 

Loss coefficient for flows through stator 

ZETAA2 

/a2 

Loss coefficient for flows through rotor 

ZETAl 

fi 

Loss coefficient for velocities through 
stator 

ZETA2 

f 2 

Loss coefficient for velocities through 
rotor without after-expansion 
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1 18 


PROGRAM PERFORM 
DIMENSION YY(4,53 
COMMON Cl»C2»C3,C 
READ 281,PO»TO» GA 
READ 282 » RPMA,DM1 
READ 283, AME,PE.,Z 
DO 118 K * 1,4 
READ 284,J,(NMAXI 
NMAX = NMAX(J) 
READ 285 (YY(J,I) 
CONTINUE 
AJ = 778. 17 
G = 32.174 
RTO = SQRTFITO) 
RFR ='WDOT*RTO/PO 
Cl=SQRTF(R/32.174 
C2=SQRTF(2.*GAM*3 
C 3=lG AM-1. )/(2.*3 
C4=SQRTF(GAM*32.1 
EXP1=(GAM-1. )/GAM 
EXP2=(GAM)/(GAM-1 
EXP3=1./(GAM-1.) 
CP = R*GAM/(AJ*{3 
PE = PE/PO 
TEIS = PE**EXP1- 
DTE = ( 1.-TEIS)*( 
TE - l.-OTE 
VE = C2 *SQRTF(DTE 
CALL CURVE (1,HS1 
CALL CURVE (2,HR? 
DO 200 J = 1,5 
DRPM = 1000 . *FLOA 
RPMA = 13000.-DR? 
PRINT 250 

251, WDOT,PO 

252, RFR,GAM 
253»C1» EXP1 

254, C2,EXP2 

255, C3,EXP3 

256, C4,RTO 
RPMA/RTO 
257 

258, RPMA,R? 

259, DM1,HS1 

260, DM2,HR? 

261, ALFAB1 

262, RETAB 1 

263, BET AB2 

264, ZET AA2 
= ALFAB1*3 


),NMAX(4),TI(4 ) 

4, GAM, EXP I,EXP2,EXP3,P0,T0,R 
M,R,WDOT . 

.DM2.HS1,HR2,ALFAB1,BETAB1,BET A82,AE2 
ETAAl,ZETAA2»AMIS 

J ), TI (J ) ) 

,1=1,NMAX) 


1 1 1 


ioo 


101 


103 


PRINT- 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
RPM = 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
ALFAB1 
BET AB 1 = BETAB1*3 
BETAB2 = 8ETAB2*3 
PRINT 265 
PRINT 266,AME,VE, 
L = 2 

DO 193 I = 1,30 
PI = 470. 

PI = Pl-FLOATF(I) 
PI = P1/P0 
CONTINUE 
PRAT=P1/PE 
P 2 = PI 

IF (P1-PE) 100,10 
CALL AFTER (AME,P 
IF (ZETAA) 116,11 
GO TO 105 
CONTINUE 
ALFA1=ALFAB1 
ZET AA=0.0 
CONTINUE 

CALL CURVE (3,AMI 
ZETA1 = ZETAS*HKS 
TRAT = P1**EXP1 


) 

2. 174*R/(GAM-1.)) 
AM*32.174»R) 

74*R) 

. ) 

AM-1. ) ) 

1.-ZETAA1) 

) 

,HKS,TI,NMAX,YY) 
»HKR,TI,NMAX,YY) 

TF(J-l) 

M 

, TO 

»R 


M 


.1415927/180. 

.1415927/180, 

.1415927/180, 

PE, TE 


* 10 . 


1,100 

RAT,ALFAB1,ALFA1,ZETAA) 
6, 103 


S,ZETAS,TI,NMAX,YY) 
+ZETAA 
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OTRAT = (1.-TRAT) 
T1 = 1.-OTRAT 
VI = C2*SQRTF(DTR 
AMV1 ■ VI/(C4+SQR 
IF (L-3) 119,120, 

120 PRINT 290 
L = 1 

119 CONTINUE 
PRINT 267 
AL = ALFA 1*180./J 
PRINT 291,AMV1,VI 
CALL TRNGL (1,V1, 
DTR = C 3*W1«W1 
TR 1 =T 1 + DTR 
PR P = (TR1/T1 ) **E 
PR1 = P1*PRP 
AMW1 = W1/(C4*SQR 
PRINT 269 
BE = BETA 1 * 180./3 
PRINT 268,AMW1,Wl 
IF(AMW1-1.0)104,1 

104 W1=W1*C0SF(BETAB1 
PRINT 271 

DTR = C 3«W1*W1 
TR1=T 1 + OTR 
PRP = (TR1/T1 )**E 
PR1 = P1 * PR P 
AMW1 = Wl/(C4*SQR 
BET A1 = BETAB1 
BE = BETA 1»180./3 
PRINT 268,AMW1,Wl 
GO TO 106 

105 CONTINUE 

CALL BEFORE (AMWl 
PRINT 270 
Wl = W1Y 

AMWl = W1/(C4*SQR 
PR1 = PR 1Y 
BET A1 = BETAB1 
BE = BETAl*180./3 
PRINT 268,AMWl,Wl 
IF (ZETAB) 114,13 

114 PRINT 288 
GO TO 198 

106 CONTINUE 
UV = U1/V1 

CALL CURVE _ 

ZETA2 = ZET AR»HKR 
PPRAT=PR1/P2 
TRATR = 1./(P PR AT 
DTE2 = (TR1-TRATR 

. TE2 = TR1-0TE2 
WE = C2*SQRTF(DTE 
AMWE = WE/(C4*SQR 
PRINT 272 
CALL CPHI (ZETAAZ 
PRINT 273,PHIC,PH 
PRINT 274,PPRAT,P 
IF (PPRAT-PPCR) 1 

107 CONTINUE 

' DELTA = ABSF(PPRA 
IF (DELTA-.0005)1 

122 PE2 = P2 
PRINT 275 
PRINT 266,AMWE,W= 
PRINT 276 
DT2=(TR1-TRATR+TR 
T2 = TR1-DT2 
W2=C2*SQRTF(DT2) 
AMW2 = W 2/(C 4 * S 0R 
BETA2=BETAB2 • 
BE2= BETA2*180./3 
PRINT 291,AMW2,W2 


*(1.-ZETA1) 
AT) 

TF(T 1 )) 

120 


.1415927 

,P1,T1,AL,ZETA1 

DM1,ALFA1,RPM,VM,VU1,WU1,Ul,Wl ,BETA1) 


XP2 

TF(T 1 ) ) 

.1415927 
,PR 1,TR1,BE 
04,105 
-BETA1) 


XP2 

TF(TD) 

.1415927 
,PR 1,TR1,BE 

,BET AB1,BETA1,P1,Wl,TR1,T1,W1Y,PR 1 Y, T 11 
TF(TIY)) 


.1415927 
,PR1,TR1,BE 
6,106 


(4,UV,ZETAR,TI,NMAX,YY) 


**EXP1) 

«TR1)»(1.-ZETAA2) 

2 ) 

TF(TE2)) 

,PPRAT,PHIC,PHICR,PPCR) 

ICR 

PC R 

07,115,108 

T-PPCR) 

15,115, 122 


, PE2 ,TE2 
1)*(1.-ZETA2) 

TF(T2 ) ) 

.1415927 

,P2,T2,BE2,ZETA2 


,ZE T AB ) 
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108 
1 10 

121 
1 15 


299 


112 

123 

109 


1 16 
1 17 


113 


CALL TRNGL (2,W2,DM2,BETA2,RPM,WM,WU2,VU2,U2,V2,ALFA2) 
PRINT 277 

AMV2 = V2/(C4*SQRTF(T2)) 

AL2= ALFA2*180./3.1415927 
PRINT 268,AMV2,V2,P2»T2,AL2, 

GO TO 113 
CONTINUE 
PRINT 286 
L = 6 

P2 = P2 + . OO'Ol 
PPRAT = PR1/P2 

CALL CPHI (ZETAA2,PPRAT,PHIC,PHICR,PPCR) 

DELTA = ABSF(PPRAT-PPCR) 

IF (DELTA-.0005)121 ,121,1 10 

PI = P2 

GO TO 111 

CONTINUE 

AMWE •= 1.00 

PE2 = P2 

PPRAT = PR1/PE2 . 

TRATR = 1 ./(PPRAT**EXP1) 

DTE2 = ( rRl-TRATR*TRl )*( 1.-ZE.TAA2) 

TE2 = TR1-DTE2 
WE = C2*SQRTF(DTE 2) 

AMWEX= WE/(C4*SQRTF(TE2) ) 

FORMAT (1OOX,F10.3) 

PRINT 275 

PRINT 266,AMWE,WE,PE2,TE2 
DT2=(TR1-TRATR*TR1 )»( 1.-ZETA2) 

T2 = TR1-DT2 
W2=C2*SQRTF(DT2) 

RET A2 =EE T AB2 

CALL TRNGL (2,W2, DM2,BETA2,RPM,WM,WU2,VU2,U2,V2,ALFA2) 
AMW2 = W2/(C4*SQRTF(T2)) 

AMV2 = V2/(C4*S0RTF(-T2)) 

PRINT 276 

BE2= BETA2*180./3.1415927 

PRINT 291,AMW2,W2,P2,T2,BE2,ZETA2 

PRINT 277 

AL2= ALFA2*180./3.1415927 
PRINT 268,AMV2,V2,P2,T2,AL2 
L = 3 
GO TO 113 
P2 = P2-.01 

IF (L — 5) 109, 123, 123- 

PRINT 290 

L = 1 

CONTINUE 

PRINT 287 

PRAT=P2/PE2 

CALL AFTER (AMWE,PRAT,BETAB2,BETA2,ZETAA) 

IF (ZETAA) 116,116,117 

PRINT 289 

GO TO 199 

CONTINUE 

PPRAT = PR1/P2 

TRATR = 1./(PPRAT**EXP1) 

DT2=(TR1-TRATR*TR1)*(1.-(ZETA2+ZETAA)) 

T2 = TR1-DT2 
W2=C2*SQRTF(DT2) 

CALL TRNGL (2,W2, DM2,BETA2,RPM,WM,WU2,VU2,U2,V2,ALFA2 ) 
AMW2 = W2/(C4«SQRTF(T2)) 

AMV2 = V2/(C4*SQRTF(T2)) 

PRINT 276 

BE 2= BETA2*180./3.1415927 

ZET = ZET A2 + ZET AA 

PRINT 291,AMW2.W2,P2,T2,BE2.ZET 

PRINT 277 

AL2= ALFA2*180./3.1415927 
PRINT 268,AMV2,V2,P2,T2,AL2 
CONTINUE 
P2T = P2*P0 
T2T = T2* TO 
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DTF = C 3* V2 »V2 
TT2 = T2+0TF 

TP2 = P2T*((TT2/r2)**EXP2) 
PRATO = 1./P2 
TT2 = TT2»TO 
OHl = T0-TT2 

DH2 = TO* (1.-{(TP2/PO)**EXPl)) 
EFF = DH1/DH2 
VU1T = VUl*RTO 
VU2T = VU2*RT 0 


HP = RTO*WDOT*l . + 14 * ( VUl T.*U 1-VU2T*U2 ) /(G* AJ) 

PRINT 273 

PRINT 279,VUlT,VJ2T,P2T,TP2,T2T, TT2 
PRINT 280,PRATO,HP,EFF 
L = L + 1 

IF (AMWE-1.00) 198,112,112 

2500F0RMAT { 1 H1 ,2 1X 19HTURBI NE PERFORMANCE/ / 1 3X3 5Fi INPUT 
IT ANTS COMPUTED /) 

251OFORMAT(13H FLOW RATE = F6.1, 17H LBM/SEC PO = F6.0, 
1T0 = F6.0. 2H R ) 

2520F0RMAT (22H REFERRED FLOW RATE * F5.3,4X8HGAMMA = 

/) 

= F9.4,21X7HEXP1 
= F9.4,21X7HEXP2 
= F9.6,21X7HEXP3 
= F9.4,21X7HRT 0 
PHYSICAL PARAMETERS 


VALUES AND CONS 


13H PSIA. 


1 F6.2 

253 FORMAT 
FORMAT 
FORMAT 
FORMAT 


F4.2»7X5HR 


254 

255 

256 


(7X5HC1 
(7X5HC2 
(7X5HC3 
(7X5HC4 


258 


FORMAT!/20X20H 
FORMAT!7H.RPM 


F6.4) 
F6.4) 
F6.4) 
F6.2) 
/) 


3/) 

IN 


= F5. 2»27H IN ROTOR BLADE HEIG 


262 

263 

264 

265 

266 

2670F0RMAT 
1T1/TO 

268 FORMAT 

269 FORMAT 
2700F0RMAT!29H 

10 TRI/TO 
271OFORMAT!15H 
1ETA1 ) 

272 FORMAT!27H 


) 

) 

/) 


=7.0,23X15HREFERRED RPM = F12 
2590F0RMAT(27H MEAN OIA AT ROTOR INLET = F5.2.27H 
1GHT = F5.3,3H IN ) 

2600F0RMATI27H MEAN OIA AT ROTOR EXIT 
1HT = F5.3.3H IN ) 

261 FORMAT!12X26H STATOR DISCHARGE ANGLE = 

FORMAT! 12X26H ROrOR INLET ANGLE 
FORMAT! 12X26H ROrOR DISCHARGE ANGLE = 

FORMAT ( 7X41HL0SS COEFFICIENT THROUGH 
FORMAT!18H STATOR EXIT PLANE//32H M(VE) 

FORMAT (F6.2.F10.3.2F8.4/) 

!2 3H STATOR 01SCHARGE PLANE//49H M!V1) 

ALFA1 ZETA1 ) 

(F6.2,F10.3,2F8.4,F8.2/) 

( 41H M(W1 ) Wl/RTO PR 1/PO TRI/TO BETA! 

IN ROrOR BEHIND ENTRY SH0CK//41H M(W1) 

BETA1 ) 

IN ROTOR ENTRY//41H M(Wl) Wl/RTO PR1/P0 


STATOR BLADE HE I 


F6.2,4H DEG 
f6.2,4H DEG 
F6.2,4H DEG 
ROTOR (FLOWS) 
VE/RTO 


= F7.4//) 
PE/PO TE/TO) 


/1/RTO 


PI /PO 


) 


Wl/RTO PR1/P 


TRI/TO 


B 


FLOW FUNCTION COMPUTATIONS) 

273 FORMAT (10X16H PHI COMPUTED = F6.5,13X.15HPHI .CRITICAL = F6.5) 
2740F0RMAT!26H PRESSURE RATIO IMPOSED = F6.3,28H CRITICAL PRESSURE RA 
1 TI 0 = F6.3) 

275 FORMAT (17H ROTOR EXIT PLANE//32H M!WE) WE/RTO PE2/P0 TE2/T0) 
2760F0PMAT (22H ROTOR DISCHARGE PLANE//49H M(W2) W2/RT0 P2/P0 T 
12/TO BETA2 ZETA2 ) 

277 FORMAT(41H M(V2) V2/RT0 P2/P0 T2/T0' ALFA2 ) 

2780F3RMAT ( /15X31H RESULTANT POWER AND EFFICIENCY // 

1 3X11HVU1!FT/SEG ),2X11HVU2(FT/SEC),3X8HP2(PS IA),2X9HP2T(PSI A),3X 

2 2HT2,5X3HT2T) 


279 FORMAT 

280 FORMAT 
1CIENCY 

281 FORMAT 

282 FORMAT 

283 FORMAT 

284 FORMAT! 

285 FORMAT 

286 FORMAT(/48H 

287 FORMAT (21H 

288 FORMAT 

289 FORMAT 

290 FORMAT 

291 FORMAT 

198 CONTINUE 

199 CONTINUE 


(F12.2.F13.2,FI 1.0,2F10.0,F7.0/) 

(3X25HCVERALL PRESSURE RATIO = F6.3,2X5HHP = F6.0,15H 
= F 4.3 / / /) • 

(5F10.0) 

(9F8•0) 

(5F10.0) 

13,16,F10. 0) 

(6F9.0)_ 

CANNOT SATISFY FLOW FUNCTION, INCREASE PI /) 
AFTEREXPANSION) 

IS NEGATIVE. ) 

IS NEGATIVE. ) 


EFFI 


ROTOR 

ROTOR 

ZETAB 

ZETAA 


( 1 9H 
( 1 9H 
( 1H1 ) 

(F6.2.F10.3,2F8.4,F8.2,F8.4/) 
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200 


ALFA81 = 
BETAB1 = 
BETAB2 = 
CONTINUE 
END 


ALFA81*180./3.1415927 
BETAB1*180./3.1415927 
BETAB2*!80./3.1415927 


“129“ 


ggjgoug?i::s 


T' '’O* rrp 


| oucpift . 

lie-, T.^,of, * v 1 


S =-( -•>) ’’ 

in. = o 


= &OC - £ 


Confute / Y 


JLl. L'. a 


F. 
c 

Pe/i 


To 


P, , 

a 1 + ^ ? 


c / 1 V-o 


S 
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i 

L: 
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SUBROUTINE AFTER 



FORTRAN NAMES 

, EQUIVALENT SYMBOLS AND DEFINITIONS 

ALFA 

COSRAT 


Undeflected flow ang^fT (blade angle) 

Cosine ratio of blade angle to flow angle 

DALFA 

Ac* 

c/g “ 0( (degrees) 

DALFAR 

A0( 

of b ' (radius) 

PRAT 

(p d /p e > 

Imposed pressure ratio 

PRAT I 


Computed pressure ratio 

VD 

v d 

Velocity at discharge plane 

VELRAT 

V V e 

Velocity ratio of flow at discharge plane 
to exit plane 

ZA 


Loss coefficient due to after-expansion 
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SUBROUTINE AFTER ( AME,PRAT.ALFAS,ALFA 1, ZA ) 

COMMON Cl»C2»C3,C4,GAM,EXP1,EXP2»EXP3,P0»T0,R 
ALFA = ABSF(ALFAS) 

EE= (2./((GAM-1.) *AME*AME)) 

DALFA = 0.0 
M « 0 

400 M = M+l 

401 OALFA = (-.1)**M-1J+DALFA 
DAIFAR = DALFA*3.1415927/180. 

COSRAT = COSF(ALFA)/COSF(ALFA-DALFAR) 

AA = GAM*AME*AME*COSF(DALFAR)-COSRAT*(GAM-1.)*AME*AME/2. 

BB = -GAM*AME*AME-1. 

CC = COSRAT*{l.+((GAM-1.)*AME*AME/2.)) 

DO = 8B*8B-4.*AA»CC 
IF (DO) 401,402 ,402 

402 VELRAT = (-BR+SQRTF(DO))/(2.*AA) 

OPRATI = (COSRAT*( 1. + {(GAM-1.)*AME*AME/2.)*(1.-VELRAT*VELRAT) )) 
1 /VELRAT 

IF (PRAT 1) 410,411,411 

410 PRINT 409,N,AME,PRAT,ALFA,DALFAft 
409 FORMAT (72X,13,4F10.5) 

411 CONTINUE , * 

ZA = 1.-(VELRAT*VELRAT)'/(1.+EE-EE*(PRAT1**EXP1) ) 

VD = VE*VELRAT 

DIFF = ABSF(PRAT1-PRAT) 

IF (DIFF-.00001) 408,408,403 

403 IF (PRAT 1-PRAT) 404,408,401 

404 M=M+1 

405 DALFA = ( - . 1 ) »* (-1 ) +D ALF A 
DALFAR = DALFA*3.1415927/180. 

COSRAT = COSF(ALFA)/COSF(ALFA-DALFAR) 

AA = GAM*AME*AME»COSF(DALFAR)-COSRAT*(GAM-1.)*AME*AME/2. 

BB = -GAM*AME*AME-1• 

CC = COSRAT*(l.+( (GAM-1.)*AME*AME/2.)) 

DD = BB*BB-4.*AA*CC 
IF(DD) 400,406,406 

406 VELRAT = (-BB+SQ^TF(OD))/(2.*AA) 

OPRAT 1 = (COSRAT*( l.+ ((GAM-1.)*AME*AME/2.)*(1.-VELRAT*VELRAT))) 
1 /VELRAT 

IF ( PRAT 1 ) 412,413,413 

412 PRINT 409,N,AME,PRAT,ALFA,DALFAR 

413 CONTINUE 

ZA = l.-(VELRAT*/ELRAT)/{1.+ EE-EE*(PRAT 1**EXP1)) 

VD = VE*VELRAT 

DIFF = ABSF(PRAT1-PRAT) 

IF (DIFF-.00001) 408,408,407 

407 IF(PRAT-PRAT 1 ) 400,408,405 

408 CONTINUE 

IF (ALFAS) 414,415,415 

414 ALFAl = ALFAS+DALFAR 
RETURN 

415 ALFAl = ALFAS - DALFAR 

END . •. 
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SUBROUTINE TRNGL 


FORTRAN NAMES. EQUIVALENT SYMBOLS AND DEFINITIONS 


A 

B 

N 

U 

VM 

VU1 

WU1 

W1 


m 


ul 


W 


ul 


W, 


Input angle of flow 

Output flow angle (If A = o( , B 
B ■ 0< .) 



and if A 



Value selects plus or minus SIGN 
Peripheral velocity of blade at mean diameter 
Meridional component of absolute velocity 
Peripheral component of absolute velocity 
Peripheral component of relative velocity 
Relative velocity 
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SUBROUTINE TRNGL (N,V1,DM,A,RPM,VM,VU1,WUl»U,Wl,B) 
COMMON Cl,C2,C3,C4,GAM,EXP1,EXP2,EXP3,PO,TO,R 
SIGN = {- 1 ) »*N 
U = 3.1415927*RPM*DM/720. 

VM = Vl#COSF(A) ' 

IF (A) 701,702,732 

701 VUl = -V1 * SINF(A) 

WUl = VUl+SIG N*U 
B = ATATMF ( WU 1 /VM) 

W1 = WUl/SINFIB) 

702 CONTINUE 

VUl = V1 •* SI NF ( A ) 

WUl = VUl+SIGN*U 
B = ATANF(WUl/VM) 

W1 = WU1/SINFIB1 
END 



-135= 


f^ch’O'PIj 



= 136- 





















SUBROUTINE BEFORE 



FORTRAN NAMES 

, EQUIVALENT SYMBOLS AND DEFINITIONS 

A 


Interim computation for velocity ratio 

ALFA 

^R 

Blade angle entering rotor 

ALFAI 


Flow angle ahead of rotor 

AMD 

M d 

Mach number of flow after shock 

B 


Interim computation for velocity ratio 

cc 


Cosine ratio of angles involved 

D 


2 

A - B, used to determine whether shock should 
be computed, or the component of the velocity 
in direction of blade be used. 

DALFA 

A<X 

Ao< = 0(' 1 - <* R 

DTR 

‘ T rr T i) 

Temperature difference between relative total 
temperature and static temperature 

FF 


Interim computation to find Mach number after 
shock 

GG 


Interim computation to find Mach number after 
shock 

PP 

p l /p l* 

Static pressure ratio through shock 

PRP 

P Rl /p l 

Pressure ratio through shock 

PR1Y 

P Rl* /P o 

Relative total pressure behind shock 

P1Y 

p l* /P o 

Static pressure behind shock 

W 

W 1 */W 1 

Velocity ratio of flow ahead of and behind 
shock 

ZETAB 

f 

' B 

Loss coefficient of flow due to entry shock 
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SUBROUTINE BEFORE (AM 1,ALFA,ALFA1, P1 , W1,TR1,T 1,W1Y,PR1Y,T1Y,ZETAB) 

COMMON C1,C2,C3,C4,GAM,EXP1,EXP2,EXP3,PO,TO,R 

OALFA = ALFA1-ALFA 

c: = COSF(ALFA1J/COSFCALFA) 

A = ( l./(GAM+l.y)*( ( l./(CC*AM'1»AM1) )+GAM*COSF(DALFA) ) 

B = -(1(2./(GAM+1.))*{(1./(AM1*AM1))+GAM)) 

D = A *A-R 
IF (D) 602,60 1., 63 1 
602 W1Y = WHCOSFIALFA1-ALFA) 

ZETAB = 0.0 

DTR = C3*WlY*WlY 

T1Y = TR1-DTR 

PRP = (TR1/T1Y)**EXP2 

PR 1Y = P1Y* PRP 

RETURN 

601 CONTINUE 

VV = A-SORTF(A*A-B) 

PP=CC*((1 .+(GAM-1 .)«AM1*AM1/2.)*(1./VV) — ((GAM-1.)* AMI *AM 1*VV/2 . ) 1 
FF = (GAM-1.)*AM1*AM1/2. 

GG = l.-VV*VV 

AMD = VV*AM1*SQRTF(1./(1.+FF*GG)) 

ZETAB = 1.-VV *VV-(2./((GAM-1•1 *AM 1*AM1))»(((PP)**E XP11-1.) 

W1Y = W1 * VV 

P1Y = P1*PP 

T1Y = TR1-C3*W1Y*W1Y 

PRP = (TR.1/T1Y)#*EXP2 

PR 1Y = P1Y*PRP 

END 
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SUBROUTINE CURVE. 


AA 

RB 

cc 

I 

J 

K 

M 

N 

NMAX 

ORD 

T1 

T2 

X 

XA 


FORTRAN NAMES, EQUIVALENT SYMBOLS AND DEFINITIONS 

\ 


a 

b 


Tl/2 


The terms of equation ORD = a + bx + cx 


Subscripts tor ordinate data points 


Counter to locate three nearest data paints 
to input abscissa 

Array selector 

Number of data points in selected array 

Ordinate computed, i,e M loss coefficient 
value desired 

Data abscissa interval 

One-half of data interval Used to select the 
three nearest points 

Abscissa in terms of three nearest points 
Input abscissa 
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SUBROUTINE CURVE!N,XA1,ORD,TI,NMAX,YY) 

DIMENSION YY(4,50),NMAX(4),TI(4),Y(501 
COMMON Cl ,C2,C3,:4,GAM,EXPl,EXP2, EXP3,P0,T0,R 
XA = XA1 
NMAX=NMAX(N) • 

T I = T I ( N ) 

DO 306 I = 1 , NM AX 
306 Y(I)= YY(N,I) 

M = XA/TI 
X = T I *FLOATF(M) 

M = M +1 ' 

IF (M-l) 304,300, 301 

304 PRINT 212 
RETURN 

300 I = M 

J = M+l 
K = M+2 
GO TO 305 

301 NM = NMAX-1 

IF ( M-NM) 302,503,303 
303 I = M-l 
J = M 
K = M + l 
GO TO 305 

302 01FR = XA-X 
T2 = TI/2. 

IF (0IFR-T2) 303,303,300 

305 CONTINUE 

A A = Y{ I ) , 

CC = 0.5«(Y(K)+Y( I )-2.*Y{J)) 

BB = Y( J ) —Y (I )-C: 

XA = (XA/TI)-FLOATF(1-1) 

ORD = AA+B8*XA+C:*XA*XA 

212 FORMAT (36H ABSCISSA WRONG IN SUBROUTINE CURVE ) 
END • • 

END 
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SULHOuTIIlE OPHI 
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SUBROUTINE CPHI 


A 

CPHI 

DAB 

DAN 

DIFF 

EN 

EXP1 

EXP4 

EXP6 

EXP7 

GAM 

PHI 

PHI 2 

PHICR 

PPCR 

RC 

RA 

RB 

ZETA 

ZETAP 


FORTRAN NAMES, EQUIVALENT SYMBOLS AND DEFINITIONS 


n 




<# 



(P/P) 


cr 


f 

fp 


Interim critical pressure ratio while computing 

Compute subroutine name 

Name for the term (p/P) exp ^ 

Name for the term (p/P) ex ^ 

Difference between given and computed loss 
coefficient 

Polytropic exponent 

Exponent #1 = ( W - 1)/ ^ 

Exponent #4 3 (1 - ^ p ) ( W ~ 1)/^ 

Exponent #6 = 2/n 

Exponent #7 = (n+l)/n 

Specific heat ratio 

Flow function computed 

Flow function, squared 

Choked flow value of flow function 

Choked flow value of pressure ratio 

Computed loss coefficient 

Name for the term (p/P) eX ^ 

Name for the term (p/P) eX ^ 

Loss coefficient 

Polytropic loss coefficient 
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SUBROUTINE CPHI (ZETA,A,PHI,PHICR,PPCR) 

COMMON Cl»C2,C3»C4,GAM»EXP1,EXP2»EXP3»P0»T0,R 
Z6TAP = -.001 
502 Z6TAP = .OOUZETAP 

EXP4 • (1.-ZETAP)•(GAM-1.l/GAM 
RA = 1 ./( A*»*FXP4) 

R 1 ? = 1«/(A**EXP1) . 

RC = (RA-RB)/(1.-RB) 

DIFF = ZETA-RC 
IF (DIFF) 501,500,502 
501 ZETAP = ZETAP-.03001 

EXP4 = (1.-ZETAP)*(GAM-1.)/GAM V 

RA * 1 •/(A**EXP4) 

RB * 1./(A**EXP1) 

RC = (RA-RB)/(l.-RB) 

DIFF * ZETA-RC 
IF (DIFF) 501,503,500 
500 EN * GAM/(l.+ZETAP*(GAM-1•)) 

EXP6 = 2./EN 
EXP7 = (EN+1.)/EN 
DAN = l./( A**EXP6.) 

DAB - 1./(A**EXP7) 

PHI 2 * 2.*(GAM/(3AM-1.))•(DAN-DAB) 

PHI = SORTF(PHI 2) 

PHICR * ((2./(EN*l.))**(l./(EN-l.)))*SQRTF(2.»GAM*(EN-l.)/ 
1 ((GAM-1.)*(EN+1.))) 

PPCR = ((EN+1.)/2.)**(EN/(EN-1.)) 

END 
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PROGRAMS TURPLOT AND TURBINE 


AKE 

AKIS 

ALFAl 

ALFA2 

BETAl 

BETA2 

DBETA 

DETA 

EFF 

ETA 

ETAA 

ETAC 

EXX 

NALFAl 

NPHIE 

PHI 

PHIR 

PSI 


FORTRAN NAMES, EQUIVALENT SYMBOLS AND DEFINITIONS 




Leaving loss coefficient 


Head coefficient, eq. 1 

0{ ^ Angle of flow exiting stator (radians), 

absolute 


0 (2 Angle of flow leaving rotor (radians), 

absolute 




1 


Flow angle entering rotor (degrees), 
relative 


A 

Y 

*7 

7 

?A 

% 

X 

*1 

h 

f 

t 


Flow angle leaving rotor (degrees), 
relative 

¥ ' A ' A 
Y ■ 7 • ?a 

Reheat factor, defined in Ref. 1 

Efficiency 

Efficiency assumed 

Efficiency computed as an interim step for 
comparison 

Expression defined by eq. 5a 

Angle of flow exiting rotor (degrees), 
absolute 

Carryover coefficient of kinetic energy into 
next stage 

Velocity coefficient in stator 

Carryover coefficient of kinetic energy from 
stator into rotor 

Velocity coefficient in rotor 
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PSIA 

RADRAT 

RBETA1 

RBETA2 

RELRAT 

RSTAR 

VELM1 

VELRAT 

WSQR1 

WSQR2 

XSQR 


7 A 

r 2 /r x 

A 

A 


V .l' U l 


< V „2 /V »1> 


w 2 
1 


w, 


Velocity coefficient in rotor computed for 
comparison 

Radius ratio through rotor 

Flow angle entering rotor (radians), relative 

Flow angle leaving rotor (radians), relative 

Ratio of relative velocities 

Degree of reaction, that fraction of the 
expansion through the stage which occurs in 
the rotor 

Meridional component of absolute velocity 
leaving stator/Peripheral velocity of rotor 
at mean diameter 

Velocity ratio of meridional components 
The relative velocity entering rotor, squared 
The relative velocity leaving rotor, squared 
Expression defined by eq. 5a 
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PROGRAM TURPLOT 

DO 50 NALFAl * 60 ,75,5 * 

do 50 I * 1,12 

RSTAR*-.05+(FLOAT F( I) )*.05 

IF(N-2)21,22,22 • '• 

21 PRINT 40 

N = N+1 • 

GO TO 23 

22 N = 1 

PRINT 47 ' 

23 CONTINUE 
NPHIE * 0 
RADRAT = 1,00 
VELRAT =1.00 
PHI = 0.95 
PHIR = 0.95 
EFF = O.C 
ETAC = .40 
PSI * .99 

ALFA1 = FLOATF(NALFA1)*3.1415927/180. 

AKIS = 2.0 
J = 1 

13 AKIS = AKIS +2.0 
J = J + 1 

170RBETA1 * ATANFtTANFlALFA!) - ((1.0)/(PHI«CO$F(ALFA 1)+SQRTF 
1 (1. - RSTAR)»S3RTF(AKI$) ) ) ) 

BETA 1 * RBETAW130./3.1415927 
12 CONTINUE • 

OXSGR = (RST AR* (1.+EFF)) + ((PHI + PHI)*(1,-RSTAR)*(PHIR-(VELRAT* 

1 VELRAT*C0SF(ALFA1)»C0SF(ALFA1}/(PSI*PSI)))) - (2.*PHIR*PHI 

2 *SINF{ALFA!)»SQRTF(1,-RSTAR)/SQRTF(AKIS))+((PHIR+RAORAT* 

3 RAORAT-1.)/AKIS) 


19 


EXX = SQRTF (XSQR ) 


RBETA2=ATANF(-PSI *EXX/(PH I*VELRAT*COSF(ALFAl)*SQRTF(1,-RSTAR) ) ) 
BETA2 = RBETA2+130./3.1415927 
20 DBETA * BETA1-BETA2 

PSIA =' 0.99 - .030223»DBETA - 4.97/(180.-OBETA) 

DPSI = ABSFIPSI -PSIA) 

IF (.0000001-DPSI) 10,11,11 

10 PSI = PSIA 
GO TO 12 

11 PSI = PSIA 

OETA = (2./SQRTFIAKIS))»((PHI*SINF(ALFAl)*SQRTF(1.-RSTAR) )- 
1 (RADRAT*RADRAT/SQRTF(AKIS))+(PSI*EXX*RAORAT)) 

OAKE = AKIS*((PHI* 3 HI*(1.-RSTAR)*COSF(ALFA 1)*COSF(ALFA 1)*VELRAT* 

1 VELRAT) + ((PSI»PSI)*XSQR)-(2.»PSI«EXX*RADRAT/SQRTF(AKIS) ) 

2 +( (RADRAT*RADI AT)/AKIS)) 

VELM1 = PHI«SQRT=(1.-RSTAR)*SQRTF(AKIS)*COSF(ALFA)) 

0ALFA2 = ATANF(((RADRAT/SQRTF(AKIS))-(PSI*EXX))/(PHI*VELRAT* 

1 COSF(ALFAl)*S3RTF(1•-RSTAR))) . 

ALFA2 = ALFA2*lS0./3.1415927 

OWSQRl = 1•♦(4KIS» PHI*PHI*(1•-RSTAR))-2«* SQRTF(AKIS)*PHI* 

1 SQRTF(1.-RSTAR)*SINF(ALFAl) 

0WSQR2 = (PSI*PSI*AKIS)*(((1.+EFF)*RSTAR)+(PHIR*PHI«PHl« 

1 (1 .-RSTAR))-(2.»PHlR*PHI*SQRTF(l.-RSTAR)*SINF(ALFAl ) )/ 

2 (SQRTF(AKIS))+(PH IR+RADRAT*RADRAT-1.)/AKIS) 

RELRAT = SQRTF(W>QR2)/SQRTF(WSQR1) 

IF (J-6) 14,14,13 
15 PRINT 47 
J = 2 

140PRI NT 41,NALFAl,NPHIE,RADRAT,VELRAT,PHI,PHIR,PSI,RSTAR,AKIS, 

1ET A,AKE,VELM1,RELRAT,BETA1,BETA2,ALFA2 
IF (ETA-ETAC) 50, 50,13 

400F0RMAT(/I18H1 Al PE R2/R1 • VM2/VM1 PHI PHIR PSI R* 

1 KIS ETA KE VM1/U1 W2/W1 BET Al BETA2 ALFA2 /) 

410F0RMAT ( 13,16,F7.1,F9.1,F8.2,F7.2,2F7.3,F8.3,F9.5,2F7.3, 

1 F8.3,2F8.2,F7.2 ) 

47 FORMAT (/) 

50 CONTINUE 
END 
END 
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/ (PHI*COSF(ALFA 1 )«SQRTF 


PROGRAM TURBINE 

00 50 N*1,4 ' ‘ 

N A L F A 1 = 80 - ( 5* N ) 

PRINT 40 
00 50 1=1,12 

RSTAR*-.05+(FL0ATF(I) )••05 
IF (1-7) 18,20,15 
20 PRINT 40 
18 CONTINUE 

J * 1 • 

K = 1 
NPHIE * 0 

ALFAl * FLOATFtNALFAl>*3.1415927/180. 

RADRAT = 1.00 
VELRAT = 1.00 
PHI = 0.95 
PHIR = 0.95 
AKI$*8.0 
EFF * 0.0 
PSI = .99 
ETAA = 0.0 

170RBETA 1 * ATANF(TANF(ALFA1) -.(( 1.0) 

1 (1. - RSTAR)*SQRTF(AKIS)))) 

BET A1 = RBETAl«130./3.141S927 

12 CONTINUE 

OXSQR * (RSTARM1.+EFF))♦((PHI*PHI)*(1.-RSTAR)*(PHIR-(VELRAT* 

1 VELRAT*COSF(AL= A1)*COSF(ALFAl)/(PSI*PSI)))-(2.*PHIR*PHI 

2 *SINF(ALFAl)*SQRTFt1,-RSTAR)/SQRTFtAKIS))+((PHIR+RADRAT* 

3 RADRAT-1.)/AKIS) ) 

EXX = SQRTF(XSQR) 

RBETA2=ATANF(-PSI *EXX/(PH I*VELRAT*COST (ALFAl)*SQRTF(1.-RSTAR))) 
BETA2 = RBETA2*130./3.1415927 
DBETA * BETA!-BETA2 

PSIA = 0.99 - .000228*DBETA - 4.97/(180.-DBETA) 

DPSI = ADSF(PS I -PSIA) 

IF (,0000001-DPSI) 10,11,11 

10 PSI = PSIA 

K = K + 1 • • 

IF (K-1000) 12,12,45 

11 PSI = PSIA 

OETA = (2./SQRTFIAKIS) )•( (PHI *^S INF (ALFAl) *SQRTF( 1.-RSTAR) )- 
1 ( RADRAT *R ADRAT/ SQRTF(AKIS) ) ■♦■ ( PS I »EXX*RADRAT ) ) 

IF (ETA-ETAA) 13,14,14 

13 DETA = ABSF(ETA-ETAA) 

IF (.0000001-DETA) 15,16,16 

14 CONTINUE 
ETAA = ETA 

DKIS = (-.1)**J . 

AKIS * AKIS+OKIS 
GO TO 17 

15 J = J+1 

IF (J-5) 14,14,43 
43 PRINT 44 
GO TO 50 
45 PRINT 46 
GO TO 50 

16 CONTINUE 

OAKE =AKIS*((PHI»>HI*(14-RSTAR)*COSF(ALFAl)*COSF(ALFAl)*VELRAT* 

1 VELRAT)♦((PSI* PSI)*XSQR)-(2.»PSI*EXX*RADRAT/SQRTF(AKIS ) ) 

2 +((RAORAT*RADRAT)/AKIS)) 

VELM1 = PHl*SQRT=(1.-RSTAR)*SQRTF(AKIS)*COSF(ALFAl) 

0ALFA2 = ATANF(((RADRAT/SGRTF(AKIS))-(PSI*EXX))/(PHI*VELRAT* 

1 COSF(ALFAl)*SQRTF(1.-RSTAR))) 

ALFA2 * ALFA2*18D./3.1415927 

0WS3R2 = (PSI*PSI*AKIS)*(((1.+EFF)*RSTAR)+(PHIR*PHI*PHI* 

1 (1.-RSTAR))-(2.*PHIR*PHI*SQRTF(1.-RSTAR)*SINF(ALFA1))/ 

2 (SORTF(AKIS) )«• (PHIR+RADRAT*RAQRAT-1. J/AKIS] 

QWSQrT; 1«MAKIS*PHI*PHt*(1.-RSTAR) i • SQRtF (AK l S )*PHl# 

1 SQRTFt1.-RSTAR)*SINF(ALFAl) 

RELRAT = SQRTFtWSQR2)/SQRTF(WSQR1) 

OPRINT 41,NALFAl,NPHIE,RADRAT,VELRAT,PHI,PHIR,PSI,RSTAR,AKIS, 

1 ETA,AKE.VELM1,RELRAT,BETA 1,BETA2,ALFA2 
L = 1 
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19 CONTINUE 

ETAR * ETA-.OOS 

24 AKIS=AKIS+.1 
K = 1 

ORB ETA 1 = ATANF(TANF(ALFA 1) - {( 1.0) / (PHI*COSF(ALFA 1)*SORTF 

1 (1. - RSTAR)*SQRTF(AKIS)))) 

BETA 1 = RBETAl*130./3.1415927 

22 CONTINUE 

OXSQR = (RSTAR*(1.+EFF)) + ((PHI*PHI)*(1.-RSTAR)*(PHI R-tVELRAT* 

• 1 V£LRAT*COSF(ALFAl)*C0SF(ALFA11/(PS I*PSI)))-(2.*PHIR*PHl . 

2 *SINF(ALFAl)*SQ3TF(1.-RSTAR)/SQRTF(AKIS))♦((PHIR+RADRAT* 

3 RAORAT-1.J/AKIS)) 

EXX = SQRTF(XSQR) 

RBETA2=ATANF(-PSI*EXX/(PH I*VELRAT*COSF(ALFAl)*SQRTF(1.-RSTAR))) 
BETA2 = RBETA2*130./3.1415927 
DBETA * BETA!-BET A2 

PSIA = 0.99 - .0D0228*DBETA - 4.97/(180.-OBETA) 

DPSI * ABSF(PS I -PSIA) 

IF (,0000001-DPSI.) 21,25,25 
21 PS I=P SI A 
K =K+ 1 

IF (K-1000) 22,22,45 

25 PS I=PSIA 

OETA * (2./SQRTF(AKIS) )*((PHI *S INF( ALFA! ) *SQRTF (1 .-RSTAR ))- 
1 (RADRAT«RADRAT/SQRTF!AKIS)) + (PS I*EXX*RADRAT)) 

IF (ETA-ETAR) 23,24,24 

23 CONTINUE 
M= 1 

55 AK I S= AKIS-.OO’l 
M=M + 1 

IF (1000-M) 51,35,35 

35 CONTINUE 

0RBETA1 = ATANF(TANF(ALFAl) - (( 1.0) / (PHI*COSF(ALFA 1)*SQRTF 
1 (1. - RSTAR)*S3RTF(AKIS)))) 

BETA1 = RBETAl*130./3.1415927 •> 

OXSQR = (RSTAR*(1.+EFF))+((PHI*PHI)*(1.-RSTAR)*{PHIR-(VELRAT* 

1 VELRAT*COSF(ALFAl)*COSF(ALFAl)/(PSI*PSI)))-(2.*?HIR*PHI 

2 *SINF(ALFAl)*SQ3TF(1.-RSTAR)/SQRTF(AKIS))+((PHIR+RADRAT* 

3 RADRAT-1.J/AKI5)) 

EXX - SQRTF{XSOR) 

RRETA2=ATANF(-PSI• EXX/(PH I*VELRAT*COSF(ALFA1)*SQRTF(1.-RSTAR)) ) 
BETA2 = RBETA2*130./3.1415927 
OBETA = BETA1-BETA2 

PSIA = 0.99 - .000228*DBETA - 4.97/(180.-DBETA) 

DPSI = ABSF(PSI -PSIA) 

IF (,0000001-DPSI) 36,38,38 

36 PS I =P SIA 
GO TO 35 

38 PS I =PSIA 

OETA = (2./SQRTF(AKIS))*((PHI*SINF(ALFAl)*SQRTF(l.-RSTAR))- 
1 (RADRAT*RADRAT/SQRTF(AKIS)) + (PS I*EXX*RADRAT)) 

IF (ETA-ETAR) 55,26,26 

260AKE =AKIS*((PHI**HI*(1 .-RSTAR)*COSF(ALFA 1)*COSF( ALFA1)*VELRAT* 

1 VELRAT ) + ( ( PSI* PSD *XSQR)-( 2. *PS I *E XX*RADRAT/SQR TF ( AKIS) ) 

2 +((RAORAT*RADRAT)/AKIS)) 

VELM1 = PHI*SQRTF (1.-RSTAR)*SQRTF(AKIS)*C0SF(ALFAl ) 

0ALFA2 = ATANF((RADRAT/SQRTF(AKIS))-(PSI*EXX))/(P3I*VELRAT* 

1 COSF(ALFA!)*SQRTF(l.-RSTAR))) . 

ALFA2 = ALFA2*180./3.1415927 

0WSQR2 = (PSI*PSI*AKIS)*(((1.+EFF)*RSTAR)+(PHIR*PHI*PHI* 

1 (1.-RSTAR))-(2.*PHIR*PHI«SQRTF(1.-RSTAR)*SINF(ALFA1))/ , 

2 (SQRTF(AKIS))f(PHIR+RADRAT*RADRAT-1.)/AKIS) 

0WSQR1 = l.+(AKIS*PHI*PHI*(l.-RSTAR))-2.*SQRTF(AKIS)*PHI* 

1 SQRTF(1.-RSTA1)*SINF( ALFAl ) 

RELRAT = SQRTFtWSQR2)/SQRTF(WSQR1) 

OPR INT 41,NALFA1,N PHIE,RAORAT,VELRAT,PHI,PHIR,PSI»RSTAR,AKIS, 

1 ETA, AKE,VELM1,RElRAT,BETA 1,BETA2,ALFA2 
L=L + 1 

I F(L-3) 19,42,42 
42 PRINT 47 
GO TO 50 
51 PRINT 52 
GO TO 50 
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400F0RMAT(/118H1 Al PE R2/R1 VM2/VM1 PHI PHIR PSI R* 

1 KIS ETA KE VM1/U1 W2/W1 BET Al 8ETA2 ALFA2 /) 

410 FORMAT (/ I 3,1 6, F7.1,F9.1,F8.2,F7.2,2F7.3,F8.3,F9.5,2F7.3, 

1 F8.3,2F8.2,F7.2 ) 

44 FORMAT (/35H FIVE ITERATIONS OF AKIS COMPLETED) 

46 FORMAT (/34H 1003 ITERATIONS OF PSI COMPLETED ) 

52 FORMAT (35H 1000 ITERATIONS .OF AXIS COMPLETED ) 

47 FORMAT (///) 

40 FORMAT (1H1) 

50 CONTINUE 
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PROGRAM FANNQ3 


FORTRAN NAMES. EQUIVALENT SYMBOLS AND DEFINITIONS 


A 

DAB 

DAN 

EE 

EN 

EXP1 

EXP2 

EXP3 

EXP4 

GAM 

PHI 

PH12 

RA 

RB 

R 


P/P 


n 


4 

<t > 2 


ZETA 

ZETAP 


1 

f. 


Pressure ratio 

Name for the term (p/P) 6X ^ 

Name for the term (p/P) ex ^ 

Reheat factor defined by ( ^ ^ )/(1 - ^ ) 

Polytropic exponent. See Appendix I 
Exponent #1 • (2/n) 

Exponent #2 » (n+l)/n 

Exponent #3 * (1 - ^ p ) ( 'fl'- l)/^ 

Exponent #4 ■ ( # - 1 )/ / 3f' 

Specific heat ratio 
Flow function 

Flow function squared. Intermediate step in finding 

Name for the term (p/P) ex *^ 

Name for the term (p/P) ex ^ 

Computed loss coefficient to be used for 
determining p . R ■ (RA-RB)/(1-RB) 

Loss coefficient ■ 1 

Polytropic loss coefficient defined in Appendix I 
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11 

12 


13 


102 


TAPE 

TAPE 


PROGRAM FANN03 
DIMENSION PHI(11) 
DO 22 I * 125,143 
XI * I 

GAM * XI/100. 
100F0RMAT (26H1F10W 
19H PRESSURE ) 

31 OOFflRMAT (30H1 2ETA 
19H PRESSURE ) 

2100F3RMAT (32H1 REHE 
19H PRESSURE ) 
WRITE OUTPUT TAPE 
WRITE OUTPUT TAPE 
WRITE OUTPUT TAPE 
DO 11 M = 1,11 
ZETA(M) * -.025+. 
FORMAT (18H RATI 
WRITE OUTPUT TAPE 
WRITE OUTPUT 
WRITE OUTPUT 
•A * 1.02 
CONTINUE 
DO 111 N * 1,11 
ZETA(N) * -.025+. 
ZETAP(N) * -.001 
ZETAP(N) « +.001+ 
EXP3 * (1 .-ZETAP( 
EXP4 * (GAM-1.)/G 
RA * 1«/(A++EXP3) 
RB * 1./(A++EXP4) 
R « (RA-RB)/(l.-R 
DIFF * ZETA(N)-R 
IF (OIFF) 101,100 
ZETAP(N) * -.0003 
EXP3 = (l.-ZETAP( 
EXP4 = (GAM-1.1/3 
RA = 1./(A++EXP3) 
RB « l./(A*»EXP4) 
R = (RA-RB)/(l.-R 
DIFF = ZETA(N)-R 
IF (OIFF) 101,100 
EN * GAM/ ( 1. ♦ ZE T<V 
EXP 1 « 2./EN 
EXP2 * (EN+1. )/EN 
DAN » 1•/(A**EXP1 
DAB * 1•/(A++EXPZ 
PHI 2 = 2.♦(GAM/(G 
PHI(N) * SQRTF(PR 
EE(N) * (ZETAP(N) 
WRITE OUTPUT TAPE 
WRITE OUTPUT TAPE 
WRITE OUTPUT TAPE 
FORMAT (F6.2.13X. 
IF (A-1.99) 17,19 
A = A+.02 
GO TO 13 

IF (5.9-A) 22,20, 
A = A+.20 
GO TO 13 
CONTINUE . 

FORMAT (8FT2.5) 
ENO 
END 


101 


100 


111 


15 

17 

19 

20 
21 
22 

400 


,ZET A(11),ZETAP(11),EE( 1 


FUNCTION FOR GAMMA » F5, 

POLYTROPIC FOR GAMMA « 

AT FACTOR *E* FOR GAMMA 

3,10,GAM 
4,310,GAM 
5,210,GAM 

025*FL0ATF(M) 

0 ZETA » 11F9.3/) 

3, 1 2 rZE.TA 

4.12, ZETA 

5.12, ZETA . 


025«FL0ATF(N) 

ZET APIN) 

N))»(GAM-1.)/GAM 
AM 


Q) 

,102 

1+ZETAPlN) 

N))*(GAM-1.)/GAM 
AM 


B) 

,100 

P(N)»(GAM-1.)) 


) 

) 

AM-1.))*(DAN-DAB) 

12 ) 

-ZETA(N))/(l.-ZETA(N)) 

3.15, A,PHI 

4.15, A,ZETAP 

5.15, A,EE 
11F9.5) 

, 19 


20 



ID 

2 // 

F5.2// 

« F5.2// 


-152- 


PROGRAM FANNQ4 


FORTRAN NAMES. EQUIVALENT SYMBOLS AND DEFINITIONS 

A Interim critical pressure ratio while computing 


CPHI 


Compute *p , subroutine name 

DAB 


Name for the term (p/P) ex ^ 

DAN 


Name for the term (p/P) 6X ^ 

DIFF 


Difference between given and computed loss 
coefficient 

EN 

n 

Polytropic exponent 

EXP1 


Exponent #1 * ( - 1)/^ 

EXP4 


Exponent #4 ■ (1 - - 1)/^ 

EXP5 


Exponent #5 ■ 2/n 

EXP6 


Exponent #6 ■ (n+l)/n 

GAM 

I? 

Specific heat ratio 

PH 

(p/p) cr 

Critical pressure ratio name for storing away 
and printing 

PHI 

<p 

Flow function assumed for comparison 

PHI 2 


Flow function, squared 

PHIC 

to 

Flow function computed in subroutine based on 
and 5 

PHICR 

A* 

Choked flow value of flow function 

PPCR 

< P /*»cr 

Choked flow value of pressure ratio 

PR 

P/P 

Pressure ratio 

R 


Computed loss coefficient 

RA 


Name for the term (p/P) ex ^ 

RB 


Name for the term (p/P) 65 ^ 
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z 

4 

Polytropic loss coefficient name for storing away 
and printing 

ZETA 

f 

Loss coefficient 

ZETAP 

fp 

Polytropic loss coefficient 
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PROGRAM FANN04 

DIMENSION ZETA(11) ,PR(11),PH(11 ) , Z ( 1 1 ) 

PRINT UOO 
DO 499 J=1,11 

499 ZETA(J)=-.025+.025+FLOATF(J) . 

PRINT 401»ZETA 

DO 490 I = 125*140 
X I = I 

GAM= XI/100. 

DO 498 K=1,11 
A = 1 .8 
P HI = 0.2 

497 CONTINUE 

CALL CPHI (GAM,ZETA(K),A,PHIC,PHICR,PPCR,ZETAP) 

IF (PHIC-PHI) 494,493,495 

495 CONTINUE 
PHI = PHIC 
A = PPCR 
GO TO 497 

494 DIFF = ABSF(PHIC-PHI) 

IF (DIFF-.0000001) 493,493,496 

496 CONTINUE 
PHI = PHIC 
A = PPCR 
GO TO 497 

4.93 CONTINUE 

PR(K) = PPCR 
PH( K) = PHICR 
UK) = ZETAP 

498 CONTINUE 

IF( 1-133) 492,491 ,492 

491 PRINT 400 
PRINT 401,ZETA 

492 PRINT 402,GAM 
PRINT 412,PH 
PRINT 403,PR 
PRINT 404,Z 

4000F3RMAT (90H1 CH0<6D FLOW VALUES OF THE FLOW FUNCTION, PRESSURE RA1 
110 AND POLYTRCPi: LOSS COEFFICIENT //) 

401 FORMAT (11X7HZETA = 11F9.3//J 

402 FORMAT (7H GAMMA* F4.2) 

412 FORMAT (18H PHI * 11F9.5/) 

403 FORMAT (18H (PO/P)C * 11F9.5/) 

404 FORMAT (18H ZETAP = 11F9.5//) 

490 CONTINUE 

END 

SUBROUTINE CPHI (GAM,ZETA,A,PHIC,PHICR,PPCR,ZETAP) 

ZETAP = -.001 

500 ZETAP=0.001+ZETA> 

EXP4*(1.-ZETAP)*(GAM-1.)/GAM 

EXP 1 = (GAM-1. J/3AM : 

RA = 1./(A**EXP4) 

RB = 1 •/(A**EXP1) 

R = (RA-RB)/(1.-RB) 

IF(ZETA-R) 501,532,500 

501 ZETAP=-.00001+ZETAP 

EXP4=(1.-ZETAP)*( GAM-1. )/GAM 
EXP 1 = (GAM-1.)/3AM 
RA = 1./(A»*EXP4) 

R3 = 1./(A**EXP1) 

R = (RA-RB)/(l.-RB) 

IF(ZETA-R) 501,532,502 

502 EN=GAM/(1.+ZETAP*(GAM-1.)) 

EXP5 * 2./EN 

EXP6 = (EN + 1.)/EN 
DAN * 1 •/(A*»EXP5) 

DAB = 1./(A»* fcXPS) 

PHI2 = 2.*(GAM/(3AM-1.))*(DAN-DAB) 

PHIC = SQRTF(PH12 ) 

OPHICR = ((2./(EN+1. ))**(1./IEN-1.)))*SQRTF((2./EXP1)*((EN-1.)/ 

1 (EN+1.))) 

PPCR * ((EN+1.)/2.)**(EN/(EN-1•)) 

END 

* 
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100 

101 


200 


1 02 

103 

104 


BB 

CC 

00 

IF 


105 


106 

300 


PROGRAM AFTER 
AME * 1.37 
NALFA = 75 

ALFA = FLOATF(NALFA)»3.1415927/180. 

GAMMA • 1.37 

EXP = (GAMMA -1.)/GAMMA 

EE = (2./((GAMMA-1.)*AME*AME)) 

PRINT 409 

PRINT 410,AME,NALFA,GAMMA 

PRINT 404- 

DO 50C L = 1,26 

PRAT = 1.02 - .02 *FLOATF(LI 

DALFA = 0.0 

M * 0 

M = M+l 

CALFA = (-.1 ) *»(M-1l + OALFA 
DALFAR * 0ALFA*3.1415927/180. 

COSRAT = COSF (ALF A J./CCSF ( ALFA-0 ALF AR) 

AA = GAMMA*AME*AME*C0SF(0ALFAR)-C0SRAT*(GAMMA-1.)•AME*AME/2. 

* -GAMMA*AME*AME-1 • 

= COSR AT *(l. + ( (GAMMA-1•)*AME*AME/2 • ) 1 
= BB*SB-4.*AA*CC 
(OC) 101,200,200 
VELRAT1 = (-BB+SQRTF(DO 1)/(2.*AA) 

OPRAT 1 = (COSRAT»( 1.♦((GAMMA-1.1*AME*AME/2.)•(1.-VELRAT 1*VELRAT1 ) ) 3 
1 /VELRAT1 

2A = 1.-(VELRAT1*VELRAT1)/(1. + EE-EE*(PRAT1«*EXP)) 

DIFF = ABSF(PRAT1-PRAT) . 

IF (DIFF-.0001 )300,300.102 
IF (PRAT 1-PRAT) 1 03,300,101 
M = M+l 

OALFA = (-. 1 )**(M-l)+OALFA 
DALFAR = DALFA*3.1415927/180. 

COSRAT = COSFULF A)/CCSF( ALFA-DALFAR) 

AA = GAMMA*AME»AME*C0SF(CALFAR)-C0SRAT*(GAMMA-1.)»AME«AME/2. 

= -GAMMA*AME*AME-1. 

= COSRAT * ( 1 . +( (GAMMA-1.)*AME«AME/2.)) 

= BB*BB-4. *AA*CC 
(OC) 100,105,105 
VELRA.T 1 = (-BB+SQRTFtOD) ) / (2. • AA) 

OPRATI = (COSRATM l. + ((GAMMA-1. )•AME*AME/2.1•(1.-VELRAT1*VELRAT 11 )) 
1 /VELRAT1 

ZA = 1.-(VELRAT1*VELRAT1)/(1.+EE-EE•(PRAT1«*EXP)) 

DIFF * ABSF(PRATl-PRAT) 

IF (DIFF-.0001) 300,300,106 
IF (PRAT-PRAT 1 ) 100.300,1 04 

PRINT 4 12,P(^AT,PRAT1,VELRAT1 ,DALFA,ZA,AMO,AMA 
GO TO 
FCRMAT 


BB 

CC 

DO 

IF 


401 
4 020F0RMA T 


1 


4 03 

404 

405 

406 
409 
4 10 

411 

412 
500 


1 


AME 

FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
MO 

FORMAT 
CONTINUE 
ENO 
ENO 


B2-4AC' 

A 


500 
(2 OH 
(60H 
) 

( 6 F1 0.6/ ) 

(46H PRAT 

(/E20.10) 

(F7.4,4F10.4) 
(25H1 MACH 
(F10.3,I5,F10.3 
(66H PRAT 

MA /) 
(7F10.4) 


IS NEGATIVE 
B 


/) 


PRAT 1 VELRAT1 


COSRAT 


OALFA 


GAMMA 


ZA /) 


ALFA 

/) 


GAMMA 


/) 


PRAT1 VELRAT1 


OALFA 


ZA 
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Appendix IV. Critical Mach Number of Flow Undergoing Polytropic 
Expansion. 


The critical pressure ratio for a flow undergoing a polytropic 
expansion, as developed in Appendix I, is 



(1,13) 


To find the maximum velocity corresponding to this pressure 
ratio it must be recalled that, for a polytropic process 





'D-l 

T - C -p '” v 


constant 


where C is a constant 


(IV,1) 


(IV,2) 


Therefore 



(IV,3) 


Using the T-s diagram of Fig. 24 as a reference, it can be seen that 








and 


JL = _e_*L 

Zjj J r-/ 


% 


/ 


-<?r 


For the critical velocity then, using (1,13) 


\4 = Zf£jn r a 

ir-i 


/-) ' 






which reduces to 


I^ - ZfZ 7 '* £ ( ”-> 


If-/ \ /n +/ ) 


The acoustic velocity is 




7' er = V** ($)'"- v f T °(^r, 




which reduces to 




Therefore 


M ‘ 

KA. 


]/ * 
Y <yi, 




Stf —/ 

r -/ 


(IV,4) 


(IV,5) 
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(IV,6) 


(IV,7) 
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Also developed in Appendix I was the relation for the polytropic 


exponent n, in terms of the polytropic loss coefficient, as 


= 


V 

£ (v-0 +/ 


(1,7) 


Making this substitution, and collecting terms 

/ 


STL-/ _ / ~~^P 

V-/ 'fpK'-V-/ 


and 
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(IV,9) 


* 


P 


0.10, M 

cr 

0.20. M 


For example, for 2^ 


1.37 and 
and for 


.931, 

.854. 





















